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Abstract 

We survey the implications for new physics of the discrepancy between the LEP measure- 
ment of Rb and its Standard Model prediction. Two broad classes of models are considered: 
(i) those in which new Zbb couplings arise at tree level, through Z or 6-quark mixing with 
new particles, and (n) those in which new scalars and fermions alter the Zbb vertex at 
one loop. We keep our analysis as general as possible in order to systematically determine 
what kinds of features can produce corrections to Rb of the right sign and magnitude. 
We are able to identify several successful mechanisms, which include most of those which 
have been recently been proposed in the literature, as well as some earlier proposals {e.g. 
supersymmetric models). By seeing how such models appear as special cases of our general 
treatment we are able to shed light on the reason for, and the robustness of, their ability 
to explain Rb- 
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1. Introduction 



The Standard Model (SM) of electroweak interactions has been tested and confirmed 
with unprecedented precision over the past few years using measurements of e'^e~ scat- 
tering at the Z resonance at LEP [1] and SLC [2]. A particularly striking example of 
the impressive SM synthesis of the data came with the discovery, at CDF and DO [3], of 
the top quark with a mass which is in excellent agreement with the value implied by the 
measurements at LEP. 

The biggest — and only statistically important — fly to be found so far in the prover- 
bial SM ointment is the experimental surplus of bottom quarks produced in Z decays, 
relative to the SM prediction. With the analysis of the 1994 data as described at last 
summer's conferences [1][2], this discrepancy has become almost a 4cr deviation between 
experiment and SM theory. The numbers are: 

i^b = Ffo/Fhad = 0.2219 ±0.0017, while i?6(SM) = 0.2156. (1) 

The SM prediction assumes a top mass of rrit = 180 GeV and the strong coupling constant 
as{Mz) = 0.123, as is obtained by optimizing the fit to the data. 

There are other measurements which differ from their SM predictions at the > 2(7 
level: Rc (2.5cr), A^^{t) (2. Oct), and the inconsistency (2.4cr) between as measured at 
LEP with that obtained from A^^^ as determined at SLC [2]. In fact, since the Rc and Rb 
measurements are correlated, and because they were announced together, some authors 
refer to this as the "i?b-i?c crisis." One of the points we wish to make in this paper is that 
there is no Rc crisis. If the i?6 discrepancy can be resolved by the addition of new physics, 
one then obtains an acceptable fit to the data. In other words, Rc, as well as A^g{T) and 
A^j^, can reasonably be viewed simply as statistical fluctuations. 

On the other hand, it is difficult to treat the measured value of Rf, as a statistical 
ffuctuation. Indeed, largely because of Rb, the data at face value now exclude the SM at 
the 98.8% confidence level. If we suppose that this disagreement is not an experimental 
artifact, then the burning question is: What Does It Mean? 

Our main intention in this paper is to survey a broad class of models to determine 
what kinds of new physics can bring theory back into agreement with experiment. Since Rb 
is the main culprit we focus on explaining both its sign and magnitude. This is nontrivial, 
but not impossible to do, given that the discrepancy is roughly the same size as, though 
in the opposite direction to, the large m^-dependent SM radiative correction. The result 
is therefore just within the reach of one-loop perturbation theory. 
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Our purpose is to survey the theoretical possibihties within a reasonably broad frame- 
work, and we therefore keep our analysis quite general, rather than focusing on individual 
models. This approach has the virtue of exhibiting features that are generic to sundry 
explanations of the Z ^ bb width, and many of the proposals of the literature emerge as 
special cases of the alternatives which we consider. 

In the end we find a number of possible explanations of the effect, each of which 
would have its own potential signature in future experiments. These divide roughly into 
two categories: those which introduce new physics into Rb at tree level, and those which 
do so starting at the one-loop level. 

The possibilities are explored in detail in the remainder of the article, which has 
the following organization. The next section discusses why i?5 is the only statistically 
significant discrepancy between theory and experiment, and summarizes the kinds of in- 
teractions to which the data points. This is followed by several sections, each of which 
examines a different class of models. Section 3 studies the tree-level possibilities, consisting 
of models in which the Z boson or the b quark mixes with a hitherto undiscovered particle. 
We find several viable models, some of which imply comparatively large modifications to 
the right-handed 6-quark neutral-current couplings. Sections 4 and 5 then consider loop 
contributions to Rb. Section 4 concerns modifications to the t-quark sector of the SM. 
Although we find that we can reduce the discrepancy in Rb to ~ 2cr, we do not regard 
this as sufficient to claim success for models of this type. Section 5 then considers the 
general form for loop-level modifications of the Zbb vertex which arise from models with 
new scalars and fermions. The general results are then applied to a number of illustrative 
examples. We are able to see why simple models, like multi-Higgs doublet and Zee-type 
models fail to reproduce the data, as well as to examine the robustness of the difficulties 
of a supersymmetric explanation of Rb- Finally, our general expressions guide us to some 
examples which do make experimentally successful predictions. Section 6 discusses some 
future experimental tests of various explanations of the Rb problem. Our conclusions are 
summarized in section 7. 



2. The Data Speaks 



Taken at face value, the current LEP/SLC data excludes the SM at the 98.8% confi- 
dence level. It is natural to ask what new physics would be required to reconcile theory 
and experiment in the event that this disagreement survives further experimental scrutiny. 
Before digging through one's theoretical repertoire for candidate models, it behooves the 
theorist first to ask which features are preferred in a successful explanation of the data. 
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An efficient way to do so is to specialize to the case where all new particles are 
heavy enough to influence Z-pole observables primarily through their lowest-dimension 
interactions in an effective lagrangian. Then the various effective couplings may be fit to 
the data, allowing a quantitative statistical comparison of which ones give the best fit. 
Although not all of the scenarios which we shall describe involve only heavy particles, 
many of them do and the conclusions we draw using an effective lagrangian often have 
a much wider applicability than one might at first assume. Applications of this type of 
analysis to earlier data [4] [5] have been recently updated to include last summer's data [6], 
and the purpose of this section is to summarize the results that were found. 

There are two main types of effective interactions which play an important role in the 
analysis of Z-resonance physics, and we pause first to enumerate briefiy what these are. 
(For more details see Ref. [4].) The first kind of interaction consists of the lowest-dimension 
deviations to the electroweak boson self-energies, and can be parameterized using the well- 
known Peskin-Takeuchi parameters S and T [7].^ The second class of interactions consists 
of nonstandard dimension-four effective neutral-current fermion couplings, which may be 
defined as follows:^ 

J^e^ = T^ZJ^^ [{g[ + Sg[) 7. + {gi + SgQ 7.] /. (2) 

In this expression g[ and g{i denote the SM couplings, which are normalized so that 
g[ = l[ — and gt^ = —Q^sl^^ where Ig and are the third component of weak 

isospin and the electric charge of the corresponding fermion, /. s^, = sin^^, denotes the 
sine of the weak mixing angle, and 7i(i,) = (1 T 75)/2. 

Fitting these effective couplings to the data leads to the following conclusions. 

• (1) What Must Be Explained: Although the measured values for several observables 
depart from SM predictions at the 2cr level and more, at the present level of experimental 
accuracy it is only the i?^ measurement which really must be theoretically explained. After 
all, some 2a fluctuations are not surprising in any sample of twenty or more independent 
measurements. (Indeed, it would be disturbing, statistically speaking, if all measurements 
agreed with theory to within la.) This observation is reflected quantitatively in the flts 
of Ref. [6], for which the minimal modification which is required to acommodate the i?5 
measurement, namely the addition of only new effective Zhh couplings, already raises the 

^ The third parameter, U, also appears but doesn't play a role in the Z-pole observables. 

2 f 

Here we introduce a slight notation change relative to Ref. [4] in that our couplings Sg[^ correspond 

to Sgl of Ref. [4]. 
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confidence level of the fit to acceptable levels (x^j„/d.o.f. = 15.5/11 as compared to 27.2/13 
for a SM fit). We therefore regard the evidence for other discrepancies with the SM, such 
as the value of i?c, as being inconclusive at present and focus instead on models which 
predict large enough values for R^. 

• (2) The Significance of Re'- Since the 1995 summer conferences have highlighted the 
nonstandard measured values for the Z branching ratio into both c and h quarks, it is 
worth making the above point more quantitatively for the particular case of the discrep- 
ancy in i?c. This was addressed in Ref. [6] by introducing effective couplings of the Z to 
both h and c quarks, and testing how much better the resulting predictions fit the obser- 
vations. Although the goodness of fit to Z-pole observables does improve somewhat (with 
X^.^/d.o.f. = 9.8/9), it does so at the expense of driving the preferred value for the strong 
coupling constant up to as{Mz) — 0.180 ± 0.035, in disagreement at the level of 2a with 
low-energy determinations, which lie in the range 0.112 ± 0.003 [8]. This change in the 
fit value for as{Mz) is driven by the experimental constraint that the total Z width not 
change with the addition of the new Zee couplings.^ Once the low-energy determinations 
of as{Mz) are also included, x^.^yd.o.f. not only drops back to the levels taken in the 
fit only to efi'ective Zbb couplings, but the best-fit prediction for R^ again moves into a 
roughly 2a discrepancy with experiment. 

It is nevertheless theoretically possible to introduce new physics to account for i?c in 
a way which does not drive up the value of the strong coupling constant. As argued on 
model-independent grounds in Ref. [6], and more recently within the context of specific 
models [9] [10], an alteration of the c-quark neutral-current couplings can be compensated 
for in the total Z width by also altering the neutral-current couplings of light quarks, such 
as the s. We put these types of models aside in the present paper, considering them to be 
insufficiently motivated by the experimental data. 

• (3) LH vs. RH Couplings: The data do not yet permit a determination of whether it is 
preferable to modify the left-handed (LH) or right-handed (RH) Zbb coupling. The mini- 
mum values for found in Ref. [6] for a fit involving either LH, RH or both couplings are, 
respectively, xi.„/d.o.f.(LH) = 17.0/12, xL„/d.o.f.(RH) = 16.1/12 or xLJd.o.f.(both) = 
15.5/11. 



Introducing effective 6-quark couplings have precisely the opposite effect — since the SM prediction 
for Vb is low and that for Fc is high relative to observations — lowering the strong coupling constant to 
as(Mz)=0.103±0.007. 
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L Figure 1: A fit 

of the Zbb couphngs Sg^ to Z-po\e data from the 1995 Summer Conferences. The four sohd hnes 
respectively denote the la, 2a, 3a, and 4(7 error ellipsoids. The SM prediction lies at the origin, (0, 0). 
This fit yields as{Mz) = 0.101 ± 0.007. 



• (4) The Size Required to Explain Rb'- The analysis of Ref. [6] also indicates the size of 
the change in the neutral-current 6-quark couplings that is required if these are to properly 
describe the data. The best fit values which are required are displayed in Figure 1, and 
are listed in Table 1. Table 1 also includes for comparison the corresponding tree- level 
SM couplings, as well as the largest SM one- loop vertex corrections (those which depend 
quadratically or logarithmically on the t-quark mass^ m^), evaluated at = 0.23. For 
making comparisons we take rrit = 180 GeV. 

As we now describe, the implications of the numbers appearing in Table 1 depend on 
the handedness (LH vs. RH) of effective new-physics Zbb couplings. 

• (4-a) LH Couplings: Table 1 shows that the required change in the LH Zbb couplings must 

^ More precisely [11], we use 6g'l = {j^)[r-\-2.9,S\nr], where r=m?/M^. 
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Coupling 


^(SM tree) 


5^(SM top loop) 


5^ (Individual Fit) 


8g{¥\t to Both) 




-0.4230 


0.0065 


-0.0067 ±0.0021 


-0.0029 ± 0.0037 


A 


0.0770 





0.034 ±0.010 


0.022 ±0.018 



Table 1 

Required Neutral-Current b-Quark Couplings: The last two columns display the size of the effective cor- 
rection to the left- and right-handed SM Zhh couplings which best fit the data. The "individual fit" is 
obtained using only one effective chiral coupling in addition to the SM parameters mt and ots [Mz)- The "fit 
to both" includes both couplings. Also shown for comparison are the SM predictions for these couplings, 
both the tree- level contribution ("SM tree"), and the dominant mt-dependent one-loop vertex correction, 
evaluated at s^=0.23 ("SM top loop"). 

be negative and comparable in magnitude to the mt-dependent loop corrections within the 
SM. The sign must be negative since the prediction for the Z ^ bb width must be increased 
relative to the SM result in order to agree with experiment. This requires 5g^ to have the 
same sign as the tree- level value for g^, which is negative. As we shall see, this sign 
limits the kinds of models which can produce the desired effect. Comparison with the SM 
loop contribution shows that the magnitude required for 6g'^ is reasonable for a one-loop 
calculation. Since the size of the mt-dependent part of the SM loop is enhanced by a factor 
of m|/M^, the required new-physics effect must be larger than a generic electroweak loop 
correction. 

• (4b) RH Couplings: Since the SM tree- level RH coupling is opposite in sign to its LH 
counterpart and is some five times smaller, the new-physics contribution required by the 
data, dg^, is positive and comparable in size to the tree-level coupling. This makes it likely 
that any new-physics explanation of the data which relies on changing g^ must arise at 
tree level, rather than through loops. 

• (5) Absence of Oblique Corrections: A final proviso is that any contribution to g'l or 
g'^ should not be accompanied by large contributions to other physical quantities. For 
example, Ref. [6] finds that the best-fit values for the oblique parameters S and T are 

5= -0.25 ±0.19 

(3) 

T= -0.12 ±0.21 
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(with a relative correlation of 0.86) even when 6g'l ^ are free to float in the fit. Since T 
often gets contributions similar in size to 5g^, these bounds can be quite restrictive. 

Notice that we need not worry about the possibility of having large cancellations 
between the new-physics contributions to the oblique parameters and Sg'l in R^. It is 
true that such a partial cancellation actually happens for in the SM, where the loop 
contributions proportional to m| in T and Sg'l exactly cancel in the limit that = |, and 
so end up being suppressed by a factor — |. We nevertheless need not consider such 
a cancellation in i?6 since the oblique parameters (especially T) almost completely cancel 
between Tb and Fhad- Quantitatively, we have [4]: 

4.57 Sgl + 0.828 Sg^ - 0.00452 S + 0.0110 
1.01 5gl + 0.183 Sg''^ - 0.00518 S + 0.0114 t) (4) 
3.56 Sgi + 0.645 + 0.00066 S - 0.0004 t) . 

We now turn to a discussion of the circumstances under which the above conditions 
may be achieved in a broad class of models. 

3. Tree-Level Effects: Mixing 

At tree level the Zbb couplings can be modified if there is mixing amongst the charge — | 
quarks, or the neutral, colourless vector bosons. Being a tree-level effect it is relatively easy 
and straightforward to analyze and compare different scenarios. Also, since mixing effects 
can be large, mixing can provide comparatively large corrections to the Z66-coupling, such 
as is needed to modify Ri, through changes to g^. Not surprisingly, a number of recent 
models [9] [10], [12]-[14] use mixing to try to resolve the Ri, (and Re) discrepancy. Our aim 
here is to survey the possibilities in a reasonably general way. We therefore postpone for 
the moment a more detailed phenomenological analysis of the various options. 

In general we imagine that all particles having the same spin, colour and electric 
charge can be related to one another through mass matrices (some of whose entries might 
be constrained to be zero in particular models due to gauge symmetries or restrictions 
on the Higgs-field representations). We denote the colour-triplet, charge Q — — |, quarks 
in the flavour basis by S", and label the corresponding mass eigenstates^ by ¥. The 

We imagine having already diagonalized the SM mass matrices so that in the absence of this non- 
standard mixing one of the B" reduces to the usual b quark, with a diagonal mass matrix with the d and 
s quarks. 



n = rr (i - 

T — -PSM / 1 

i had — J- had I ~ 
SO Rb = Rb^ (l - 
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mass-eigenstate quarks, 6*, are obtained from the by performing independent unitary- 
rotations, (ul^n^ , amongst the left- and right-handed fields. The b quark that has been 
observed in experiments is the lightest of the mass eigenstates, b = b^, and all others are 
necessarily much heavier than this state. 

Similar considerations also apply for colourless, electrically-neutral spin-one particles. 
In this case we imagine the weak eigenstates, Z^, to be related to the mass eigenstates, Z^, 
by an orthogonal matrix, Al™"*. We take the physical Z, whose properties are measured in 
such exquisite detail at LEP and SLC, to be the lightest of the mass eigenstates: = Z^. 

Assuming that all of the and Z"^ (except for the lightest ones, the familiar b and Z 
particles) are too heavy to be directly produced at Z-resonance energies, we find that the 
fiavour-diagonal effective neutral-current couplings relevant for i?^ are 



2 



u 



al 



M 



wl 



(5) 



where the neutral-current couplings are taken to be diagonal in the flavour 3°^ basis. ^ 

This expression becomes reasonably simple in the common situation for which only two 

S\ fb 



particles are involved in the mixing. In this case we may write = ^ ^/ J ' ~ v ^ / 

= ^ ^/ ^ 7 ^-i^d take , I^r and A4 to be two-by-two rotation matrices parameterized 
by the mixing angles 9i^, On and 9z- In this case eq. (5) reduces to 



9l,r 



(9. 



^L,R 



+ ) 



'l,rSl,r 



+ 



l,r'^l,r 



(9 



L,R'^L,R 



(6) 



where s^, denotes sin^j^, etc. Increasing i?5 requires increasing the combination (g^)^ + 
(g^)'^- To see how this works we now specialize to more specific alternatives. 



Eq. (5) describes the most relevant effects for the Rb problem, namely the mixing of Z and b with 
new states. However, in general other indirect effects are also present, such as, for example, a shift in Mz 
due to the mixing with the Z'. For a detailed analysis of the simultaneous effects of mixing with a Z' and 
new fermions, see Ref. [15]. 
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3.1) Z Mixing 



First consider the case where two gauge bosons mix. Then eq. (6) reduces to 

9ln=i9z)L,nCz + {g^>)r.,nSz, (7) 

where {gz)L,R is the SM couphng in the absence of Z mixing, and {gz')L,R is the 6-quark 
couphng to the new field Z'^ (which might itself be generated through 6-quark mixing). It 
is clear that so long as the Z'bb coupling is nonzero, then it is always possible to choose 
the angle dz to ensure that the total effective coupling is greater than the SM one, {gz)L,R- 
This is because the magnitude of any function of the form f{9z) = Acz+B Sz is maximized 
by the angle tan^^ = B/A, for which \f\^^^ = \A/cz\ > \A\. 

The model-building challenge is to ensure that the same type of modifications do 
not appear in an unacceptable way in the effective Z couplings to other fermions, or in 
too large an Mz shift due to the mixing. This can be ensured using appropriate choices 
for the transformation properties of the fields under the new gauge symmetry, and suffi- 
ciently small Z-Z' mixing angles. Models along these lines have been recently discussed in 
Refs. [9], [16]. 

3.2) b-Quark Mixing 

The second natural choice to consider is pure 6-quark mixing, with no new neutral 
gauge bosons. We consider only the simple case of 2 x 2 mixings, since with only one 
new B' quark mixing with the SM bottom quark, eq. (6) simplifies considerably. As we 
will discuss below, we believe this to be sufficient to elucidate most of the features of the 
possible 6-mixing solutions to the problem. 

Let us first establish some notation. We denote the weak SU{2) representations of 
the SM Bl,r and of the -B^^^ as Rl,r and -R'j^^^, respectively, where R = {I,h). The 
SM S-quark assignments are R^ = (|, — |) and Rn = (0,0). By definition, a B' quark 
must have electric charge Q = —1/3, but may in principle have arbitrary weak isospin 

In terms of the eigenvalues Ig^. and of the weak-isospin generator Ig acting on B'^ 
and 5^, the combination of couplings which controls becomes 

oc {gtr + {gl? = (-f + f + ''^^-) ' + (f + '^^^«) ' • 
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In order to increase using this expression, 9^^ and must be such as to make more 
negative, more positive, or both. Two ways to ensure this are to choose 



There are also two other alternatives, involving large mixing angles or large B' representa- 
tions: /^^ > 0, with sl{r^^+\) > 1 - 28^3 ~ 0.85, and /^^ < 0, with 4 |/^^| > 2^/3 ~ 
0.15. Note that, in the presence of LH mixing, the CKM elements Vqb {q = u,c,t) get 
rescaled as Vqb Cr^Vqb, thus leading to a decrease in rates for processes in which the b 
quark couples to a VF. Therefore charged-current data can in principle put constraints on 
large LH mixing. For example, future measurements of the various t-quark decays at the 
Tcvatron will allow the extraction of Vtb in a model-independent way, thus providing a lower 
limit on Cl. At present, however, when the assumption of three-generation unitarity is re- 
laxed (as is implicit in our cases) the current measurement of BR{t — > Wb) / BR{t — > Wq) 
implies only the very weak limit \Vtb\ > 0.022 (at 95% c.l.) [17]. Hence, to date there are 
still no strong constraints on large LH mixing solutions. Regarding the RH mixings, as 
discussed below there is no corresponding way to derive constraints on c^, and so large 
solutions are always possible. 

We proceed now to classify the models in which the SM bottom quark mixes with 
other new Q = —1/3 fermions. Although there are endless possibilities for the kind of 
exotic quark one could consider, the number of possibilities can be drastically reduced, 
and a complete classification becomes possible, after the following two assumptions are 
made: 

• (i): There are no new Higgs-boson representations beyond doublets and singlets. 

• (ii): The usual S-quark mixes with a single B' , producing the mass eigenstates b and 
b' . This constrains the mass matrix to be 2 x 2: 



We will examine all of the alternatives consistent with these assumptions, both of 
which we believe to be well-motivated, and indeed not very restrictive. The resulting mod- 
els include the "standard" exotic fermion scenarios [18] (vector singlets, vector doublets, 
mirror fermions), as well as a number of others. 



or 




(9) 




(10) 
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Let us first discuss assumption (i). From Table 1 and eq. (8) one sees that the mixing 
angles must be at least as large as 10% to explain i^^, implying that the ofF-diagonal entries 
in the mass matrix eq. (10) which give rise to the mixing are of order SL^nM22 ^ O{10) 
GeV. If these entries are generated by Higgs fields in higher than doublet representations, 
such large VEVs would badly undermine the agreement between theory and experiment 
for the My^^/Mz mass ratio J 

According to assumption (i)^ the permitted Higgs representations are Ru = (|, ±|) 
and (0,0). It is then possible to specify which representations -R^^^^ allow the B' to mix 
with the B quark of the SM: 

• (1): Since the B' should be relatively heavy, we require that M22 7^ 0. Then the 
restriction (i) on the possible Higgs representations implies that 

1^1-41 = 0,^; (11) 



and 



l4.-^LI=0,i. (12) 



• (2): To have b-b' mixing, at least one of the off-diagonal entries, M12 or M21, must be 
nonzero. These terms arise respectively from the gauge-invariant products (S> Rl ® -R^ 
and Rh ® Rr ® -R^ so that -R^^^j^ must transform as the conjugate of the tensor product 



Rh ® Rr{l)'- 



or 



i?l = i?«(8)i?H = (0,0), Q,±^^ 



(13) 



i?; = i?„®i?, = (0,0), 0,-0 ,(1,-1), (1,0). (14) 

Thus the only possible representations for the B' are those with = , | , 1 and = 
0, i, 1 , |, subject to the restrictions (11)-(14). 



The contribution of these relatively large non-standard VEVs cannot be effectively compensated by 
new loop-effects. On the other hand, beyond Higgs doublets, the next case of a Higgs multiplet preserving 
the tree-level ratio is that of I3h=3, Yh='2. We do not consider such possibilities, which would also require 
the mixed B' to belong to similarly high-dimensional representations. We also neglect alternative scenarios 
invoking, for example, more Higgs triplets and cancellations between different VEVs, since these suffer 
from severe fine-tuning problems. 
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As for assumption (ii), it is of course possible that several species of B' quarks mix 
with the B, giving rise to an N x N mass matrix, but it seems reasonable to study the 
allowed types of mixing one at a time. After doing so it is easy to extend the analysis to 
the combined effects of simultaneous mixing with multiple B' quarks. Thus (ii) appears 
to be a rather mild assumption. 

There is one sense in which (ii) might appear to restrict the class of phenomena we 
look at in a qualitative way: it is possible to obtain mixing between the B and a B' 
in one of the higher representations we have excluded by "bootstrapping", that is, by 
intermediate mixing with a B[ in one of the allowed representations. The idea is that, if 
the SM B mixes with such a B[, but in turn the latter mixes with a B2 of larger isospin, 
this would effectively induce a B-B'2 mixing, which is not considered here. However, since 
mass entries directly coupling B to B2 are forbidden by assumption (i), the resulting B B2 
mixing will in general be proportional to the B~B[ mixing, implying that these additional 
effects are subleading, i.e. of higher order in the mixing angles. This means that if the 
dominant B B[ mixing effects are insufficient to account for the measured value of 
adding more B' quarks with larger isospin will not qualitatively change this situation. 

There is, however, a loophole to this argument. If the mass matrix has some symmetry 
which gives rise to a special "texture," then it is possible to have large mixing angles 
and thus evade the suppression due to products of small mixing angles alluded to above. 
Indeed, we have constructed several examples of 3 x 3 quasi-degenerate matrices with 
three and four texture zeros, for which the B-B2 mixing is not suppressed and, due to the 
degeneracy, can be maximally large. For example, let us choose B'l in a vector doublet with 
l3L,R = +1/2 and B2 in a vector triplet with Isl,r = +1- Because of our assumption of no 
Higgs triplets, direct B mixing with such a B2 is forbidden, and M13 = M^i = M12 — 0. 
It is easy to check that for a generic values of the nonvanishing mass matrix elements, the 
induced s]^^^ mixings are indeed subleading with respect to s]^ p.. However, if we instead 
suppose that all the nonzero elements are equal to some large mass /U, then there are two 
nonzero eigenvalues mf,/ ~ n and nib" ~ 2// while the B — B2 mixing angles s^^ ~ y^l/3 
and s]^ ~ •y/3/8 are unsuppressed relative to s^^^.^ Although it may be unnatural to 
have near-equality of the mass entries generated by singlet and doublet Higgs VEVs, as 
is needed in this case and in most of the other examples we found, it is still possible that 
some interesting solutions could be constructed along these lines. 

Apart from some special cases analogous to the one outlined above, we can therefore 
conclude that neither does assumption (ii) seriously limit the generality of our results. 

Q 

A small perturbation of the order of a few GeV can be added to some of the nonzero mass entries to 
lift the degeneracy and give a nonzero value for TO5. 
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Table 2 

Models and Charge Assignments 
All the possible models for B-B' mixing allowed by the assumptions that (i) here are no new Higgs 
representations beyond singlets and doublets, and (ii) only mixing with a single B' is considered. The 
presence of LH or RH mixings which can affect the b neutral current couplings is indicated under 'Mixing'. 
Subleading mixings, quadratically suppressed, are given in parenthesis. Equivalent models, for the purposes 
of i?b, are indicated by a prime (') in the 'Model' column, while models satisfying eq. (9) and which can 
account for the deviations in Ri with small mixing angles, are labeled by an asterisk Large RH mixing 
solutions are labeled by a double asterisk ^**\ while models 7, 7' and 7" allow for a solution with large 
LH mixing. 
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We can now enumerate all the possibilities allowed by assumptions (i) and (ii). 



With the permitted values of and listed above, and the requirement that at least 
one of the two conditions (13)-(14) is satisfied, there are 19 possibilities, listed in Table 
2. Although not all of them are anomaly-free, the anomalies can always be canceled by 
adding other exotic fermions which have no effect on R},. Since only the values of 1^^^ and 
/g^ are important for the b neutral current couplings, for our purpose models with the 
same /g^ ^ assignments are equivalent, regardless of ^ or differences in the mass matrix 
or mixing pattern. Altogether there are 12 inequivalent possibilities. Equivalent models 
are indicated by a prime (') in the 'Model' column in Table 2. 

Due to gauge invariance and to the restriction (i) on the Higgs sector, in several cases 
one of the off-diagonal entries M12 or M21 in eq. (10) vanishes, leading to a hierarchy 
between the LH and the RH mixing angles. If the b' is much heavier than the b, M12 = 
yields Sj, ~ M21/M22, while the LH mixing is suppressed by M^^. If on the other hand 
M21 = 0, then the suppression for is quadratic, leaving as the dominant mixing 
angle. For these cases, the subdominant mixings are shown in parentheses in the 'Mixing' 
column in Table 2. Notice that while models 2 and 6 allow for a large right-handed mixing 
angle solution of the Rb anomaly, the "equivalent" models 2' and 6' do not, precisely 
because of such a suppression. 

Six choices satisfy one of the two conditions in eq. (9), and hence can solve the Rf, 
problem using small mixing angles. They are labeled by an asterisk in Table 2. Three 
of these models (10,11,12) satisfy the first condition for solutions using small LH mixings. 
Since for all these cases I^^ < 0, a large RH mixing could alternatively yield a solution but 
because s^j is always suppressed with respect to s^, this latter possibility is theoretically 
disfavored. The other three choices (models 3,8,9) satisfy the second condition for solutions 
using small RH mixing. It is noteworthy that in all six models the relevant mixing needed 
to explain Ri, is automatically the dominant one, while the other, which would exacerbate 
the problem, is quadratically suppressed and hence negligible in the large mf,/ limit. There 
are two choices (models 5, 6) for which /g^ < and there is only RH mixing, and one (model 
2) for which /g^ < and is unsuppressed with respect to s^. These three cases allow 
for solutions with large RH mixings, and are labeled by a double asterisk ^**\ Finally, a 
solution with large LH mixing is possible (models 7, 7' and 7") in which /g^ = +1/2, and 
Igj^ = implies no RH mixing effects. 

In the light of Table 2 we now discuss in more detail the most popular models, as well 
as some other more exotic possibilities. 
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ure 2: The experimentally allowed mixing angles for a mirror family. The thick line covers the entire 
area of values for Sj^ and which are needed to agree with the experimental value for R^, to the 2(7 level 
or better. The thin line represents the one-parameter family of mixing angles which reproduce the SM 
prediction. Notice that the small-mixing solution, which passes through Sj, = Sjj = 0, is ruled out since 
= implies that any LH mixing will reduce and thus increases the discrepancy with experiment. 

• Vector singlet: Vector fermions by definition have identical left- and right-handed gauge 
quantum numbers. A vector singlet (model 1) is one for which — I'^ = 0. Inspection of 
eq. (8) shows that mixing with such a vector-singlet quark always acts to reduce Rb.^ 

^ A Q=+2/3 vector singlet can however be used to reduce Rc [10] [12] [14], provided that steps are 
taken, as suggested in Section 2 above, to avoid the resulting preference for an unacceptably large value 
for asiMz))- 
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• Mirror family: A mirror family (model 2) is a fourth family but with the chiralities 
of the representations interchanged. Because Ig^ vanishes, LH mixing acts to reduce the 
magnitude of g'l, and so tends to make the prediction for Rb worse than in the SM. For 
sufficiently large RH mixing angles, however, this tendency may be reversed. As was 
discussed immediately below eq. (9), since I'^j^ is negative a comparatively large mixing 
angle of > 1/3 is needed to sufficiently increase Such a large RH mixing angle 
is phenomenologically permitted by all off-resonance determinations of g'^ [19]. In fact, 
the 6-quark production cross section and asymmetry, as measured in the 7-Z interference 
region [21] [22], cannot distinguish between the two values s| = and 48^/3, which yield 
exactly the same rates. Hence this kind of model can solve the Ri, problem, though 
perhaps not in the most aesthetically pleasing way. As is shown in Fig. 2, the allowed 
range of mixing angles is limited to a narrow strip in the si — s\ plane. 

• Fourth family: A fourth family (model 4) cannot resolve Rh via tree-level effects because 
the new B' quark has the same isospin assignments as the SM h quark, and so they do not 
mix in the neutral current. Two other possibilities (models 4' ^^nd 4") yield the same 
^'iL,R assignments as the fourth family model, and are similarly unsuccessful in explaining 
Rh since they do not modify the h quark neutral current couplings. 

• Vector doublets: There are two possibilities which permit a Q = — | quark to transform 
as a weak isodoublet, and in both cases mixing with the SM h is allowed. They can be 
labeled by the different hyper charge value using the usual convention Q = I^, + Y . 

With the straightforward choice — 1^^^ = — 1/2 (model 5), we have Y[ = = 1/6. 
This type of model is discussed in Ref. [13], where the isopartncr of the B' is a top-like 
quark T' having charge -|-|. Since these are the same charge assignments as for the standard 
LH 6-quark, this leads to no mixing in the neutral current amongst the LH fields, and 
therefore only the right-handed mixing angle s^ is relevant for Ri,. Since /g^ is negative a 
comparatively large mixing angle of s| > 1/3 is needed to sufficiently increase Rb, in much 
the same way as we found for the mirror-family scenario discussed above. The required 
mixing angle that gives the experimental value, Rb = 0.2219 ± 0.0017, is 

4 = 0.367i°:°}i (15) 

The other way to fit a Q = —1/3 quark into a vector doublet corresponds to I^^ = 
I^^ = +1/2 (model 8) and so Y' = -5/6 [10]. The partner of the B' in the doublet is then 

The current 90% c.l. upper bound s^<0.010 [20] holds in the small mixing angle region s^<Cl/3. 
These models have the further difficulty that, except in certain corners of parameter space [23], they 

produce too large a contribution to the oblique parameters, S and T, to be consistent with the data. 
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an exotic quark, R, having Q = —4/3. Here Ig^ has the wrong sign for satisfying eq. (9) 
and so mixing decreases the magnitude of g^. On the other hand, I^j^ has the right sign 
to increase g^. Whether this type of model can work therefore depends on which of the 
two competing effects in i?b wins. It is easy to see that in this model the M21 entry in the 
B-B' mass matrix eq. (10) vanishes, which as discussed above results in a suppression of 
quadratic in the large mass, but only a linear suppression for s^. Hence, becomes 
negligible in the large mi,' limit, leaving as the dominant mixing angle in Ri,. The 
mixing angle which reproduces the experimental value for i?5 then is 

si = 0.059iH}i (16) 

However, in order to account for such a large value of the mixing angle in a natural way, 
the b' cannot be much heavier than ~ 100 GeV. 

Similarly to the Y' = —5/6 vector doublet case, models 3 and 9 also provide a solution 
through RH mixings. In model 5, the subdominant competing effect of is further 
suppressed by a smaller /g^, while in model 9 the effect of is enhanced by /g^ = +1, 
and hence a mixing angle a factor of 4 smaller that in (16) is sufficient to explain i?^. 

• Vector triplets: There are three possibilities for placing a vector B' quark in an isotriplet 
representation: I^^ = I'^^ = —1, 0, +1. The last does not allow for b mixing, if only Higgs 
doublets and singlets are present, and for our purposes, I^^ = I^j^ = (model 1') is 
equivalent to the vector singlet case already discussed. Only the assignment I^j^ = = 
— 1 (model 10) allows for a resolution of the Rb problem, and it was proposed in Ref. [12]. If 
B' is the lowest-isospin member of the triplet there is an exotic quark of charge Q = +5/3 
in the model. Again in the limit of large b' mass one combination of mixing angles (in this 
case Sn) is negligible, due to the vanishing of M12 in eq. (10). As a result, s^ plays the 
main role in Ri,. Agreement with experiment requires 

= 0.0127 ±0.0034. (17) 

Since the resulting change to g'l is so small, such a slight mixing angle would have escaped 
detection in all other experiments to date. 

Similarly to this case, models 11 and 12 also provide a solution through LH mixings. 
In model 11 the unwanted effects of are further suppressed, while for model 12 a LH 
mixing somewhat smaller than in (17) is sufficient to explain the data. 

Our analysis of tree level effects shows that both Z-mixing and 6-mixing can resolve 
the Rb discrepancy. 6-quark mixing solutions satisfying the two assumptions that (i) there 
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are no new Higgs representations beyond singlets and doublets, and (ii) only mixing with 
a single B' is relevant, have been completely classified. The list of the exotic new B' 
quarks with the right electroweak quantum numbers is given in Table 2. Solutions with 
small and Sl mixing angles are possible when the 5^ is the member with highest I^j^^ 
in an isodoublet or isotriplet, or when B'^ is the member with lowest 1^^^ in an isotriplet 
or isoquartet. In all these cases, new quarks with exotic electric charges are also present. 
Some other possible solutions correspond to I^^ < and are due to mixing amongst the 
RH 6-quarks involving rather large mixing angles, while for I^,^ = +1/2 we find another 
solution requiring even larger LH mixing. It is intriguing that such large mixing angles 
are consistent with all other 6-quark phenomenology. We have not attempted to classify 
models in which mixing with new states with very large values of I^lr can arise as a result 
of bootstrapping through some intermediate B' mixing. Under special circumstances, they 
could allow for additional solutions. 

For some of the models considered, the contributions to the oblique parameters could 
be problematic, yielding additional constraints. However, for the particular class of vec- 
torlike models (which includes two of the small mixing angle solutions) loop effects are 
sufficiently small to remain acceptable. This is because, unlike the top quark which be- 
longs to a chiral multiplet, vectorlike heavy b' quarks tend to decouple in the limit that their 
masses get large. Introducing mixing with other fermions does produce nonzero oblique 
corrections, but these remain small enough to have evaded detection. Exceptions to this 
statement are models involving a large number of new fields, like entire new generations, 
since these tend to accumulate large contributions to S and T. 

4. One-Loop Effects: t-Quark Mixing 

We now turn to the modifications to the Zbb couplings which can arise at one loop. 
Recall that this option can only explain Ri, if the LH 6-quark coupling, (7^, receives a 
negative correction comparable in size to the SM mt-dependent contributions. As was 
argued in section 2, it is the LH coupling we are interested in because a loop-level change 
in is too small to fix the discrepancy between the SM and experiment. 

The fact that the Rb problem could be explained if the mt-dependent one-loop contri- 
butions of the SM were absent naturally leads to the idea that perhaps the t-quark couples 
differently to the 6-quark than is supposed in the SM. If the t quark mixes significantly 
with a new t' quark one might be able to significantly reduce the relevant contributions 
below their SM values. In this section we show that it is at best possible to reduce the 

12 

Vectorlike models have the additional advantage of being automatically anomaly free. 
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discrepancy to ~ 2cr in models of this type, and so they cannot claim to completely explain 
the Rf, data. 

Our survey of t-quark mixing is organized as follows. We first describe the framework 
of models within which we systematically search, and we identify all of the possible exotic 
t-quark quantum numbers which can potentially work. This study is carried out much in 
the spirit of the analysis of b mixing presented in section 3. We then describe the possible t' 
loop contributions to the neutral-current b couplings. Since this calculation is very similar 
to computing the mt-dependent effects within the SM, we briefly review the latter. Besides 
providing a useful check on our final expressions, we find that the SM calculation also has 
several lessons for the more general t-quark mixing models. 

4.1) Enumerating the Models 

In this section we identify a broad class of models in which the SM top quark mixes 
with other exotic top-like fermions. As in the previous section concerning 6-quark mixing, 
we denote the electroweak eigenstates by capitals, T'*, and the mass eigenstates by lower- 
case letters, t*. To avoid confusion, quantities which specifically refer to the b sector will 
be labeled with the superscript ^ . By definition, a T' quark must have electric charge 
Q = 2/3, but may in principle have arbitrary weak isospin i?^^^ = (-f^^j^, I3l,r)- Following 
closely the discussion in the previous section, we make three assumptions which allow for 
a drastic simplification in the analysis, without much loss of generality: 

• (i): First, the usual T-quark is only allowed to mix with a single T' quark at a time, 
producing the mass eigenstates t and t'. 

• (ii): Second, for the Higgs-boson representations, we assume only one doublet and 
singlets. Additional doublets would complicate the analysis of the radiative corrections in 
a model-dependent way due to the extra diagrams involving charged Higgs bosons. 

• (in): Finally, certain T'-quark representations also contain new B' quarks. We denote 
the B'j^ and B'^ as 'exotic' whenever they have non-standard weak isospin assignments, that 
is, 1^1 ^ — i or /g^ 7^ 0. As we have already discussed, for exotic B' quarks b-b' mixing 
will modify the b neutral-current couplings at tree level, overwhelming the loop-suppressed 
t-t' mixing effects in Rf,. We therefore carry out our analysis under the requirement that 
any b-b' mixing affecting the b neutral-current couplings be absent. 

Our purpose is now to examine all of the alternatives which can arise subject to these 
three assumptions. According to (z), the T-T' mass matrices we consider are 2x2, and 
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can be written in the general form 



Due to our restriction (ii) on the Higgs sector, certain elements of this mass matrix are 
nonzero only for particular values of the T' weak isospin. Moreover, whenever belongs 
to a multiplet which also contains a Q = —1/3 quark, the and M12 entries of 
the B-B' and T-T' mass matrices are the same. In those cases in which the B' quark is 
exotic, assumption {iii) then forces us to set M12 = 0. In contrast, the M21 entries are 
unrelated - for example, the choice M21 = is always possible even if M21 7^ for the T 
and T' quarks. 

In order to select those representations, R'^ j^ which can mix with the SM T quark, 
we require the following conditions to be satisfied: 

• (1): In order to ensure a large mass for the t', we require M22 7^ 0. Analogously to (11) 
and (12), this implies 

1^1-41 = 0,^; (19) 

and 

\iL-iL\ = ^\ (20) 

• (2): To ensure a non- vanishing t-t' mixing we require at least one of the two ofi-diagonal 
entries, M12 or M21, to be non- vanishing. This translates into the following conditions on 
R'^ and R'^\ 

R', = R^®R^ = (0, 0) , ±^ , (21) 

or 

i?; = R,®R, = (0, 0) , Q, +0 , (1, 0) , (1, +1). (22) 

• (3): Whenever R'^ contains a Q = — 1/3 quark, and either B'^ or 5^ have non-standard 
isospin assignments, we require M12 = 0. This ensures that at tree level the neutral current 
h couplings are identical to those of the SM. Clearly, in the cases in which the particular 
R'^ representation implies a vanishing M21 element, imposing the condition M12 = 
completely removes all t-t' mixing. 
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We now may enumerate all the possibilities. From eqs. (19)-(22), it is apparent that as 
in the B' case the only allowed representations must have = , | , 1 and /{^ = , | , 1 , |. 

Consider first = 1 or |. In this case, from eq. (21), M21 = 0. Thus, we need 
M12 7^ if there is to be any t-t' mixing. The four possibilities for i?^ are shown in 
eq. (22). Of these, i?^ = (0,0) is not allowed since eq. (19) is not satisfied. In addition, 
R'^ = (|,^) and (1, 0) both contain exotic B' quarks {1^%^ = — | or —1) and so M12 is forced 
to vanish, leading to no t-t' mixing. This leaves i?^ = (1, 1) as a possibility, since the 
is not exotic = 0). If we choose R'^ such that I^^ = — |, then both B'^ and are 
SM-like, and b-b' mixing is not prohibited since it does not affect the b neutral current 
couplings. Thus, the combination R'^ = (|,|), i?^ = (1, 1) is allowed. 

Next consider = or |. Here, regardless of the value of I^^, M21 can be nonzero. 
Thus any i?^ representation which satisfies eqs. (19) and (20) is permitted. It is straight- 
forward to show that there are 11 possibilities. 

The list of the allowed values of Ig^ and which under our assumptions lead to 
t-t' mixing is shown in Table 3. There are twelve possible combinations, including fourth- 
generation fermions, vector singlets, vector doublets, and mirror fermions. Not all of these 
possibilities are anomaly-free, but as already noted one could always cancel anomalies by 
adding other exotic fermions which give no additional effects in i?^. 

It is useful to group the twelve possibilities into three different classes, according to 
the particular constraints on the form of the T-T' mass matrix in eq. (18). 

The first two entries in Table 3, which we have assigned to group A, correspond to the 
special case in which the -5^,^ and ^ have the same third component of weak isospin, 
hence leaving the h neutral current unaffected by mixing. Because both and appear 
in the same multiplets with T'^ and T^, two elements of the S-quark and T-quark mass 
matrices are equal: 

M12 = Mf2 , M22 = M2^2 • (23) 

As we will see, this condition is important since it implies a relation between the mixings 
and the rrit, mf mass eigenvalues. Although outside the subject of this paper, it is note- 
worthy that for these models the simultaneous presence of both b-b' and t-t' RH mixing 
generates new effects in the charged currents: right-handed Wtb charged currents get in- 
duced, proportional to the product of the T and B quark mixings s^s^. Compared to the 
modifications in the neutral currents and in the LH charged currents, these effects are of 
higher order in the mixing angles [18] [19] and, most importantly, they can only change the 
RH b coupling. But as noted above, is far too small to account for the measured R^ 
value using loop effects of this kind. Therefore the mixing-induced RH currents allowed in 
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Table 3 

Models and Charge Assignments 
Values of the weak isospin of and which, under the only restrictions of singlet and doublet Higgs 
representations, lead to nonzero t-t' neutral current mixing. The 'Model' column, labels the more familiar 
possibilities for the T' quarks: Vector Singlets, Mirror Fermions, Fourth Family and Vector Doublets. The 
other models are more exotic. 

models Ai and A2 are ineffective for fixing the Rh discrepancy, and will not be considered 
in the remainder of this paper. 

For the models in group B, the condition 

M12 = (24) 

holds. In the four cases corresponding to R'^^ = (1, 0) (models Bi, B3) and R'^^ = (1/2, 1/2) 
(models B2, Bq), an exotic 5^ quark is present in the same multiplet. Hence M12 has 
to be set to zero in order to forbid the unwanted tree-level b mixing effects. In the other 
three cases belonging to group B, T'j^ corresponds to the lowest component of non-trivial 
multiplets: i?'^ = (1,-1) (model B4) and i?'^ = (1/2,-1/2) (models S5, B7). For these 
values of 13^, M12 = is automatically ensured, due to our restriction to Higgs singlets 
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or doublets. Furthermore, these representations do not contain a 5^ quark, and no 5^ 
quark appears in the corresponding i?^. There is therefore no b-b' mixing. 

We should also remark that in model no 5^-quark appears in R'^. However, a B'^ is 
needed as the helicity partner of the present in i?^ = (1, 0). Because of our restriction 
on the allowed Higgs representations, 5^ must belong to R'^ = (1,0) or i?^ = (1/2, 1/2) , 
which in turn contain a new ^ . While the first choice corresponds to a type of 
mixing which we have already excluded from our analysis, the second choice is allowed and 
corresponds to model Bi. Following assumption (z), even in this case we neglect possible 
mixings of type Bi , when analysing ^3. 

Finally, the remaining three models constitute group C, corresponding to R'^ = (0, 0). 
In this group, is an isosinglet, as is the SM Tr, implying that only LH t-t' mixing is 
relevant. For C2 and C3, R'j^ does not contain a 5^, while for Ci the B'^ is not exotic. 
Hence in all the three cases the b neutral-current couplings are unchanged relative to the 
SM, and we need not worry about tree- level 6-mixing effects. 

4.2) t-Quark Loops Within the Standard Model 

Before examining the effect of t-t' mixing on the radiative correction to Zbb, we first 
review the SM computation. We follow the notation and calculation of Bernabeu, Pich and 
Santamaria [11] (BPS). The corrections are due to the 10 diagrams of Fig. 3. All diagrams 
are calculated in 't Hooft-Feynman gauge, and we neglect the 6-quark mass as well as the 
difference \ Vtb\'^ — 1. 

Due to the neglect of the 6-quark mass, and due to the LH character of the charged- 
current couplings, the t-quark contribution to the Zbb vertex correction preserves helicity. 
Following BPS we write the helicity-preserving part of the Z ^ bb scattering amplitude 
as 

r= - (^) 6(pi,Ai)r'^6(p2,A2)e^(9, A), (25) 

with 

r'^ = r^ + 5r'^, sr'' =^r7Li{s,r) . (26) 

ZTT 

where ST^^ represents the loop-induced correction to the Zbb vertex. I{s, r) is a dimension- 
less and Lorentz-in variant form factor which depends, a priori, on the three independent 
ratios: r = rn^/M^, s = M^/M^ and q'^/M^. For applications at the Z resonance only 
two of these are independent due to the mass-shell condition = M|. Moreover, for an 
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on-shell Z, non-resonant box-diagram contributions to e'^e — > bb are unimportant, and 
I{s,r) can be treated as an effectively gauge-invariant quantity. 




The contributions due to the t-quark may be isolated from other radiative corrections 
by keeping only the r-dependent part of /(s,r). BPS therefore define the difference 



F{s,r) = l{s,r)-I{s,0) . (27) 
25 



Given this function, the m^-dependence of the width Z — > 66 is obtained using 
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(28) 



In this last equation V.P.(s,r) denotes the mt-dependent contributions which enter F;, 
through the loop corrections to the gauge-boson vacuum polarizations. 

The function F'^'^'{s, r) is straightforward to compute. Although the resulting expres- 
sions are somewhat obscure, the special case s = reveals some interesting features which 
are also present in our new-physics calculations, and so we show the s = limit explicitly 
here. For s = 0, an evaluation of the graphs of Fig. 3 gives the following expressions: 
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The picture becomes much simpler after summing the diagrams to obtain the total SM 
contribution: 
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There are two points of interest in this sum. First, it is ultraviolet finite since all of the 
divergences oc l/(n — 4) have cancelled. This is required on general grounds since there 
can be no r-dependent divergences in /^^(s,r), and so these must cancel in F^^(s,r). 
A similar cancellation also occurs when new physics is included, provided that it respects 
the SUi^{2) X gauge symmetry and that the complete set of new contributions is 

carefully included. 

The second interesting feature of eq. (38) lies in its dependence on the weak mix- 
ing angle, Sw Each of the contributions listed in eqs. (29) through (35) has the form 
F* = (x* + J/*s^)/s^; however all of the terms involving have cancelled in the sum, 
eq. (38). This very general result also applies to all of the new-physics models we consider 
in subsequent sections. As will be proved in Section 5, the cancellation is guaranteed by 
electromagnetic gauge invariance, because the terms subleading in are proportional to 
the electromagnetic &-quark vertex at = 0, which must vanish. This gives a powerful 
check on all of our calculations. 

Rather than using complete expressions for F{s,r), we find it more instructive to 
quote our results in the limit r ^ 1, where powers of 1/r and s/r may be neglected. We 
do the same for the ratio of masses of other new particles to when these arise in later 
sections. Besides permitting compact formulae, this approximation also gives numerically 
accurate expressions for most of the models' parameter range, as is already true for the 
SM, even though r in this case is only ~ 4. In the large-r limit F^'^{s,r) becomes 



where the ellipsis denotes terms which are finite as r — > oo. Several points are noteworthy 
in this expression. 

• 1: The s-dependent term appearing in eq. (39) is numerically very small, changing the 
coefficient of Inr from 3 to 2.88. This type of s-dependence is of even less interest when 
we consider new physics, since our goal is then to examine whether the new physics can 
explain the discrepancy between theory and experiment in Ri,. That is, we want to see 
if the radiative corrections can have the right sign and magnitude to change Fi, by the 
correct amount. For these purposes, so long as the inclusion of q'^-dependent terms only 
changes the numerical analysis by factors <^ 25% (as opposed to changing its overall sign) 
they may be neglected. 

• 2: The above-mentioned cancellation of the terms proportional to when s = no 
longer occurs once the s-dependence is included. This is as expected since the electromag- 
netic Ward identity only enforces the cancellation at = 0, corresponding to s = in 




(39) 
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the present case. Notice that the leading term, proportional to r, is s- independent, and 
because of the cancellation it is completely attributable to graph (2a) of Fig. 3. All of the 
other graphs cancel in the leading term. Due to its intrinsic relation with the cancellation 
of the s^-dependent terms, the fact that only one graph is responsible for the leading 
contribution to Sg^ still holds once new physics is included. This will prove useful for 
identifying which features of a given model control the overall sign of the new contribution 
to 5gl 

• 3: Since the large-r limit corresponds to particle masses (in this case mt) that are large 
compared to M,^ and Mz, this is the limit where the effective-lagrangian analysis described 
in Section 2 directly applies. Then the function F can be interpreted as the effective Zbb 
coupling generated when the heavy particle is integrated out. Quantitatively, Sg^ is related 
to F by 

S9i = Q F. (40) 

• 4: The vacuum polarization contributions to Fb of eq. (28) have a similar interpretation in 
the heavy-particle limit. In this case the removal of the heavy particles can generate oblique 
parameters, which also contribute to F^. In the heavy-particle limit eq. (28) therefore 
reduces to the first of eqs. (4). 

4.3) Sg^ in the t-Quark Mixing Models 

We may now compute how mixing in the top-quark sector can affect the loop con- 
tributions to the process Z bb. As in the SM analysis, we set mi, = 0. In addition, 
following the discussion in the previous subsection, we neglect the s-dependence in all our 
expressions. We also ignore all vacuum-polarization effects, knowing that they essentially 
cancel in Ri,. Finally, in the CKM matrix, we set \Vid\ = \Vis\ = where i = t,t'. Thus, 
the charged-current couplings of interest to us are described by a 2 x 2 mixing matrix, 
just as in the neutral-current sector. In the absence of t-t' mixing this condition implies 
\Vtb\ = 1- 

For t-t' mixing, independent of the weak isospin of the T' , we write 




where c^, = cos^j;,, etc.. The matrices Ul^h are analogous to the b-b' mixing matrices 
defined in eq. (5) in our tree-level analysis of b mixing. 
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In the presence of t-t' mixing, the diagonal neutral-current couplings are modified: 

din = E i^tnf ^ 9'^:n + ~9U , (42) 

a=T,T' 

where i = t,t', and g^^J^ are the SM couplings defined in eq. (37). The new terms 
explicitly read 

9i = (iL - , ~9i = IL^I , (43) 

~9i = {lL-^^cl, ~9i = lLcl- (44) 

In addition, whenever the ^ has nonstandard isospin assignments, l!^^ ^ 1/2 or ^ 0, 
flavour-changing neutral-current (FCNC) couplings are also induced : 

9^,n= E 9l,nUT:J^^n^g^,n. (45) 
o=T,T' 

where i,j = t, t', and i ^ j. Here, 

9i'' = Q - 4.) , ~g'^' = /L^^Ch . (46) 

Eq. (41) determines the effective t and t' neutral-current couplings (eqs. (42)-(46)). 
However, the charged-current couplings depend on the matrix V = UjU^. Hence we need 
to consider also b mixing, since, as discussed in Sec. 4.1, in those cases in which the B' 
quark is not exotic {I^^ = —1/2, Ig^ = 0), we have no reason to require = I {i.e. no 
b-b' mixing). We then define the 2x2 charged current mixing matrix 

V = UjU^; Vtb = c,cf + s,sf, Vt'b = s,cf-c,sf, (47) 

which trivially satisfies the orthogonality conditions = V*^ V = /. In the absence of 
b-b' mixing, clearly Vtb Cl, ^t'b s^. We also note that, by assumption, whenever 
V ^ Ul v^e necessarily have I^^ = +1/2 (so that Ig^ = —1/2) in order to guarantee that 
the B'j^ is not exotic. From eqs. (43), (44) and (46), this implies that g*^ = = g*/ = 0, 
that is, the mixing effects on the LH t and t' neutral-current couplings vanish. 
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The Feynman rules of relevance for computing the Zbb vertex loop corrections in the 
presence of a mixing in the top-quark sector can now be easily written down: 



\/2M^ (48) 

zuu : W^'^ 7. + g'f^ In) + {g'^ 1. + -g'A in)] , 



Ztt : — 7p 



gi'' 1l + ~gn In 



where (f) are the unphysical charged scalars, and ti = t, t' . The vertices listed in eq. (48) 
reduce to the SM Feynman rules in the limit of no mixing. 

As pointed out at the end of subsection 4.1, in some groups of models equalities can be 
found between some elements of the T-T' and B-B' mass matrices. These have important 
consequences. In particular, once expressed in terms of the physical masses and mixing 
angles, the equalities of eq. (23) (which hold in the models of group A) can be written 

[W.MdiagWt]^, = [Wf M|i,gW|t]^^ = {K),^m,, cl; (a = 1,2) , (49) 

where Mjiag = diag [m^, mf], and we have used -^.fiag = <^i2"7-b/ (recall that we take 
777-6 = 0). Multiplying now on the left by (Z^f*^)]^^ and summing over a we obtain 

[yt Mdiag Ul\ ^2 = mt Vtb Sn + mt' Vt'b c« = . (50) 

For the models in group B, the vanishing of M12 implies no b mixing. Then V = Ul, 
and eq. (49) still holds in the limit Vtb Cli Vt'b s^. For the models in group C 
no particular relation between masses and mixing angles can be derived. For example, it 
is clear that in the 4*'* family model Ci, eq. (50) does not hold. However, for all these 
models I^j^ = 0. Hence, noting that all the couplings in eqs. (43), (44) and (46) are 
proportional to /g^, and defining r' = m^,/M^, squaring eq. (50) yields a relation which 
holds for all models in Table 3: 

~gi r = Vl, ~gi r' = - Vtb Vt'b~g^ • (51) 
This relation is used extensively in the calculation which follows. 
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b 




Figure 4: The additional Feyn- 
man diagrams which are required for models in which the t quark mixes with an exotic, heavy t' quark. 

How do we generalize the SM radiative correction to include t-t' mixing? First note 
that for each of the diagrams in Fig. 3, there is also a diagram in which all the t-quarks 
are replaced by t'-quarks. Second, there are 2 new diagrams (Fig. 4) due to the FCNC 
coupling of the Z to the t and t' . So to generalize the SM result to the case of mixing, 
three things have to be done: 

• (i): multiply eqs. (29)-(35) by V^^^ for the t contribution and for t' (with r — > r' ), 

• (ii): replace g\ j^ by the modified couplings in eq. (42), adding eqs. (43) and (44) respec- 
tively for t and t', 

• (in): include diagrams 3(a) and 3(b) (Fig. 4) corresponding to the FCNC couplings 
(eqs. (45)-(46)). 

A glance at the Feynman rules in eq. (48) shows that in the first step (z), a correction 
proportional to g2^^ , and independent of the Qi^^j^ couplings, is generated. This correction 
is common to all models in Table 3 - it appears even in the case in which the t NC couplings 
are not affected (4*'* family). In contrast, steps {ii) and {Hi) generate corrections which 
differ for different models. It is useful to recast them into two types, one proportional to 
the LH neutral current couplings (oc Vn,Vji,gL)i and the other proportional to the RH 
neutral current couplings (oc VibVjbgR)- The LH and RH corrections vanish respectively 
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for /gj^ = +1/2 and 1'^^^ = 0, when the corresponding neutral-current couphngs are not 
affected by the mixing. 

In the presence of mixing, the correction due to the diagrams of Fig. 3 involving 
internal t-quarks becomes 



2(d) 

E 

i=l{a) 



F^^{r) + F{~g*r) 



(52) 



where F^'^{r) is given by eq. (38) and 



8s2 



2r - 5 - 2r + 4 
r - 1 ^ (r - 1)2 



Inr 



(53) 



The third step (in) gives rise to a new contribution 



(54) 



Evaluating diagrams 3(a) and 3(b) (Fig. 4) we find 



/2 2 



Inr' Inr 

r — 1 r — 1 



r' — r 



-^"'V^l A + 1 + - 



■Inr' 



r' — 1 r — 1 



-Inr 



-Inr' Inr 



r' — 1 r — 1 



r' — r 



-Inr' Inr 



r' — 1 r — 1 



(55) 



• (56) 



Putting all the contributions together, for the general case we find 



3(6) 

i=l(a) j=l,2 



2 

ti6 



F«-(r,) + F(^l:',«, r,) + V^b Vt^;, F(^";, r, r') . (57) 



where tj = t, t' and r^ = r, r'. We note that due to eq. (51) all the divergent terms 
proportional to A cancel in the sum. Now, the correction bg\ = ^Xcorr to the SM 
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result can be explicitly extracted from eq. (57) by means of the relation V|j, = 1 — V^^,j,. 
Moreover, as anticipated it is possible to divide the various contributions to Xcorr into 
three different pieces: a universal correction, a correction due to LH mixing only, and a 
correction due to the RH mixing. Hence we write 

y —TP TpSM yUniv i yLH , yRH ( KQ\ 

^corr — ^ ^ — ^corr ' ^corr ' ^corr ' W°/ 

where 

Vl^, [F'^ {r') - F'"" {r)] , (59) 

VlF{glr) + Vl,F{giy) + VtbVt'bF{g'^\Ty) , (60) 

VlF{gi,r) + Vl,F{giy) + Vt,Vt>,F{g'^,Ty) . (61) 



V 

corr 



U mv 

CO 



corr 



yRH 

corr 



Using the explicit expressions for ^^'^^ and as given in eqs. (43), (44) and 
(46) above, together with relation (51) for the RH piece, these read 



yUniv _ ■\;2 fcorr(„ „l\ 
^corr — '^t'b Jl ^ ) 

X-^„ = (1 - 2/L) V,6 Vvb s,c, fr^ir, r') 



(62) 
(63) 
(64) 



with 



1 fr'(r'-6) , r'(3r' + 2),_^__, r (r - 6) r(3r + 2) 



+ ■ 



84 { r' -1 (r' - 1)2 



-In/ - 



r — 1 (r — ly 



Inr 



C2 1 i; \ rlnr +— — -r + -Inr 
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(r — 1) ir' — r) 



r — 1 {r — ly 
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(65) 



(66) 
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H — TTTi — mr 



r' - 1 (r' - ly 
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r' — r 



— mr Inr 
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, ; — -Inr' Inr 

r' — r\r' — 1 r — 1 



(67) 
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Note that a value of Vtb different from unity can be easily accounted for by using the 
unitary condition I Vt6 1 + \Vt'b\ = \Vtb\ = 1 - l^tsl + \Vtd\ in eqs. (62)-(67). 

As we have already pointed out, because of our requirement of no B-B' mixing when 
the B' is exotic, only when I'^^ = +1/2 can we have ^ Vtb, ^ Vfb- However, in this 
case X^^^ vanishes. Hence, without loss of generality, we can set the LH neutral current 
mixing equal to the charged current mixing in X^^^., obtaining 

X^orr = (1 - Vtb Vlb fr'^ir, r'), (68) 

(r,rO=^{-(r + r')+^ln^} • (69) 

From eqs. (62), (64) and (68) we see that there are only two independent mixing parameters 
relevant for the complete analysis of our problem: the LH matrix element Vtb and the RH 
mixing Sj^. Furthermore, note that as r' r, all the corrections in eqs. (65), (67) and 
(69) vanish, independent of the mixing angles. This comes about because of a GIM-like 
mechanism for all the pieces which do not depend on 1^^^. The /3^-dependent contribution 
from the RH fermions coupling to the Z vanishes in the limit r' ^ r as a consequence of 
eq. (50). 

In the limit r, r' » 1, for the functions /f°^'~(r, r') we obtain 

(70) 
(71) 

^K7)4(^v)}-(-) 

Let us now consider the numerical values of these corrections in more detail. Using 
rrit = 180 GeV, = 80 GeV, and — 0.23, eq. (38) gives a SM radiative correction of 

The question is whether it is possible to cancel this correction, thus eliminating the Rb 
problem, by choosing particular values of mf and the mixing angles. For various values of 
mt'i the value of Xcorr (eq. (58)) is shown in Table 4. 



/r>.0-5^{r'-r + 31n(^)}, 
/-(nr')-3^{-(r + r') + |l?^ln(^)}, 
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Table 4 

Dependence of the t-t' Mixing Results on m^i : This table indicates the dependence on the mass of the t' 
quark of the corrections to due to t-t' mixing, with the t mass fixed at 180 GeV. 

We see that even for mf > mt^ it is possible to choose I'^^^ I^j^, and the LH and 
RH mixing angles such that the correction is negative. So the discrepancy in Rf, between 
theory and experiment can indeed be reduced via t-t' mixing. 

Referring to the models listed in Table 4, the optimal choice for the weak isospin of 
the T' is I^^ = — 1/2 and = — 1, regardless of the value of mf. Furthermore, maximal 
RH mixing, s| ~ 1, is also preferred. However, even with these choices, it is evidently 
impossible to completely remove the Rb problem. From the above table, the best we can 
do is to take mf = 75 GeV and Vf,^ = = 0.6, in which case the total correction is 
^corr = —3.68. This leaves a 1.5cr discrepancy in R^, which would put it in the category of 
the other marginal disagreements between experiment and the SM. However, such a light 
t' quark has other phenomenological problems. In particular, CDF has put a lower limit 
of 91 GeV on charge 2/3 quarks which decay primarily to Wb [24]. Unless one adds other 
new physics to evade this constraint, the lightest t' allowed is about mf ~ 100 GeV. In 
this case, maximal LH mixing (V^^,^ — si ^ 1) gives the largest effect: Xcorr = —2.7. The 
predicted value of i?b is then still some 2cr below the measured number. 

Another possibility is that the charge 2/3 quark observed by CDF is in fact the t', 
while the real t-quark is much lighter, say mt ~ 100 GeV. Assuming small t-t' mixing. 
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and that the t' is the hghtest member of the new multiplet, the t' will then decay to Wh^ 
as observed by CDF, but the SM radiative correction will be reduced. This situation is 
essentially identical to that discussed above, in which the LH t-t' mixing is maximal, and 
mt' ~ 100 GeV: the SM value of Rh will still differ from the experimental measurement by 
about 2cr. The only way for such a scenario to work is if mt < M^. However, new physics 
is then once again required to evade the constraint from Ref. [24] . 

For all the possibilities of this section our conclusion is therefore the same: it is not 
possible to completely explain i?^ through t-t' mixing. The best we can do is reduce 
the discrepancy between theory and experiment to about 2(7, which might turn out to be 
sufficient, depending on future measurements. 

5. One-Loop Effects: Other Models 

Another way to change at the one-loop level is to introduce exotic new particles 
that couple to both the Z and the h quark. One-loop graphs involving such particles 
can then modify the Zhh vertex as measured at LEP and SLC. Recall once more the 
conclusion from Section 2: agreement with experiment requires the LH 6-quark coupling, 
to get a negative correction comparable in size to the SM m^-dependent contributions 
since loop-level changes to are too small to be detectable. 

In this section we first exhibit the general one-loop correction due to exotic new scalar 
and spin-half particles, with the goal of identifying the features responsible for the overall 
sign and magnitude of the result. We then use this general result to investigate a number 
of more specific cases. 

The answer is qualitatively different depending on whether or not the new scalars and 
fermions can mix, and thus have off-diagonal couplings to the Z boson. We therefore treat 
these two alternatives separately. The simplest case is when all Z couplings are diagonal, 
so that the one- loop results depend only upon two masses, those of the fermion and the 
scalar in the loop. Then the correction to the Zhh vertex is given by a very simple analytic 
formula, which enables us to easily explain why a number of models in this category give 
the 'wrong' sign, reducing Ft rather than increasing it. 

More generally however, the new particles in the loops have couplings to the Z which 
are diagonal only in the flavour basis but not the mass eigenstate basis, so the expressions 
become significantly more complicated. This occurs in supersymmetric extensions of the 
standard model, for example. After proposing several sample models which can resolve 
the i?fe problem, we use our results to identify which features of supersymmetric models 
are instrumental in so doing. 
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5.1) Diagonal Couplings to the Z: General Results 

We now present formulae for the correction to the Zbb vertex due to a foop involving 
generic scalar and spin- half particles. In this section we make the simplifying assumption 
that all of the Z-boson couplings are flavour diagonal. This condition is relaxed in later 
sections where the completely general expression is derived. The resulting formulae make 
it possible to see at a glance whether a given model gives the right sign for alleviating the 
discrepancy between experiment and the SM prediction for Rf,. 

The one-loop diagrams contributing to the decay Z ^ bb can be grouped according to 
whether the loop attaches to the b quark {i.e. the vertex correction and self-energy graphs 
of Fig. 5) or whether the loop appears as part of the gauge boson vacuum polarization 
(Fig. 6). For the types of models we consider these two classes of graphs are separately 
gauge invariant and finite, and so they can be understood separately. This is particularly 
clear in the limit that the particles within the loop are heavy compared to M^, since then 
the vacuum polarization graphs represent the contribution of the oblique parameters, S 
and T, while the self-energy and vertex-correction graphs describe loop-induced shifts to 
the 6-quark neutral current couplings, Sg'^ j^. 

Furthermore, although we must ensure that the oblique parameters do not become 
larger than the bound of eq. (3), eq. (4) shows that they largely cancel in the ratio i?^. 
We therefore restrict our attention in this section to the diagrams of Fig. 5 by themselves. 
The sum of the contributions of Fig. 5 is also finite as a result of the Ward identity which 
was alluded to in Section 3. This Ward identity relates the vertex-part graphs of Fig. 5a,b 
to the self-energy graphs of Fig. 5c, d. Since this cancellation is an important check of our 
results, let us explain how it comes about. 

We first consider an unbroken ^7(1) gauge boson with a tree- level coupling of gb to the 
the 6-quark. This gives rise to the familiar Ward identity from quantum electrodynamics: 
for external fermions with four-momenta p and p', 

ip-prr^=g,s{S-\p)-S-\p')), (74) 

where F'^ is the one-particle-irreducible vertex part and Sp{p) is the fermion propagator. 
If we denote the vertex-part contributions (Fig. 5a, b) to the effective vertex at zero mo- 
mentum transfer by 5gb, and the self-energy-induced wave function renormalization of the 
b quark by Zi,, then at one loop the Ward identity (74) reduces to gb{l + Zb){^ ~ = 
{gb + ^gh){^-i>'), or 

5gh - gbZb = 0. (75) 
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This last equation is the more general context for the cancellation which we found in 
Section 3; it states that the self-energy graphs (Fig. 5c, d) must precisely cancel the vertex 
part (Fig. 5a,b) in the limit of zero momentum transfer. Another way of understanding 
eq. (75) is to imagine computing the effective 6-photon vertex due to integrating out a 
heavy particle. Eq. (75) is the condition that the two effective operators h^h and h^b have 
the right relative normalization to be grouped into the gauge-covariant derivative: hph. 

But for the external Z boson, the Ward identity only applies to those parts of the 
diagrams which are insensitive to the fact that the U{1) symmetry is now broken. These 
include the 1/ (n — 4) poles from dimensional regularization, and also the contributions 
to the h neutral-current coupling proportional to s^, since the latter arise only through 
mixing from the couplings of the photon. 

We now return to the diagrams of Fig. 5. The first step is to establish the Feynman 
rules for the various vertices which appear. Since we care only about the LH neutral-current 
couplings, it suffices to consider couplings of the new particles to h^. 

>Cscalar = V (j) f^Lb + h.C. (76) 

and we write the Z coupling to / and (f) as 

hMii^ + dln)f + igs ct>%cf>] . (77) 

The couplings, g = {gl,gi,gs}, are normalized so that g = I3 — Qs^ for all fields, f^j^ 
and 0"*. 

In the examples which follow, the field / can represent either an ordinary spinor (e.g., 
t) or a conjugate spinor {e.g., f^). This difference must be kept in mind when inferring the 
corresponding charge assignments for the neutral-current couplings of the /. For example, 
the left-handed top quark has I^l = +| > so ^{ = | — and I^^ = 0, so g^ = — |s^. If 
the internal fermion were a top anizquark, however, we would instead have = — | + 1 
and g[ = +|s^. The latter couplings follow from the former using the transformation of 
the neutral current under charge conjugation: 7^1^71, <-> —'jf^'ju- 

We quote the results for evaluating the graphs of Fig. 5 in the limit where Mz (and 
of course m^) are negligible compared to m/ and M,^, since they are quite simple and 
illuminating in this approximation. It will be shown that the additional corrections due to 
the nonzero mass of the Z boson are typically less than 10% of this leading contribution. 
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(a) 



(b) 



(c) 



(d) 




Figure 5 Figure 6 

Figure 5: The one-loop vertex correction and self-energy contributions to the Zbb vertex due to fermion- 
scalar loops. 

Figure 6: The one-loop contributions to the Zbb vertex due to the gauge-boson vacuum polarizations. 



We find that 



=W^J2ri, \y f higf - g{)T{r) + {-gl + gl + gs) (a^ - T{t\ 



(78) 



where J^{r) and J^{r) are functions of the mass ratio r = mj/M^, 



(r - 1)2 
r 

[r - 1)2 



r — 1 — In r 
r — 1 — r In r 



(79) 
(80) 



denotes the divergent combination A<^ = ^j32+7 + ln(M|/47r/i2) +|, and Uc is a colour 
factor that depends on the SUc{S) quantum numbers of the fields and /. For example, 
TT-c = 1 if (/> ~ 1 or / ~ 1 (colour singlets); tt-c = 2 if / ~ 3 and ~ 3 or 6; tt-c = ^ if / ~ 3 
and ~ 8. 

The cancellation of divergences we expected on general grounds is now evident in 
the present example, because electroweak gauge invariance of the scalar interaction (76) 
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implies that the neutral-current couplings are related by 



9s +91 



9L-0- 



(81) 



This forces the term proportional to to vanish in eq. (78). As advertised the remaining 
term is both ultraviolet finite and independent of s^, which cancels in the combination 
9i - 9r- 

We are left with the compact expression 



S9r 



167r2 



(82) 



f't> 



Interestingly, it depends only on the axial-vector coupling of the internal fermion to the 
gauge boson and not on the vector coupling. The function of the masses J^{r) is 
positive and monotonically increasing, with T{r) ~ r as r — > and .^-'(oo) = 1, as can be 
seen in Fig. 7. 




5 10 15 20 Figure 7: Prom top to 

bottom, the functions J-ijn^/M^), Fji(r) = i^s(r) and Fj^{r^ which appear in the loop contribution 
to the left-handed Zhh vertex, sections 5.1 and 5.3. 



It is straightforward to generalize eq. (82) to include the effect of the nonzero Z boson 
mass. Expanding to first order in M|, one obtains an additional correction to the effective 
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vertex, 




967r2 \^m2yyo \ a;(M|/m2 - 1) + 1 (a;(M|/m2 - 1) + 1)2 j " 

(83) 

To see that this is typically an unimportant correction, consider the limit in which the 
scalar and fermion masses are equal, r = 1. Then the total correction (82) +(83) is 



(84) 



Although the M| correction can be significant if g[ = gi, the total correction would then 
be too small to explain the Ri, discrepancy, and would thus be irrelevant. 

5.2) Why Many Models Don't Work 

What is important for applications is the relative sign between the tree and one-loop 
contributions of eq. (82). In order to increase Rb so as to agree with the experimental 
observation, one needs for them both to have the same sign, and so Sg^ oc {g[ — g^) < 
in eq. (82). Thus an internal fermion with the quantum numbers of the 6-quark has 
gi—gn = —\ and would increase Rb- Conversely, a fermion like the t-quark has gt—gli = +| 
and so causes a decrease. Moreover, because the combination (^^ — g^^ is invariant under 
charge conjugation, the same statements hold true for the antiparticles: a h running in the 
loop would increase Ri, whereas a t would decrease it. 

It thus becomes quite easy to understand which models with diagonal couplings to 
the Z boson can improve the prediction for R^. Multi-Higgs-doublet models have a hard 
time explaining an R^ excess because typically it is the top quark that makes the dominant 
contribution to the loop diagram, since it has the largest Yukawa coupling, y^^ ~ 1, and 
the largest mass, to which the function T is very sensitive. However for very large tan/? 
(the ratio of the two Higgs VEV's), the Yukawa coupling of the t quark to the charged 
Higgs can be made small and that of the h quark can be made large, as in Ref. [25] . Fig. 7 
shows that, in fact, one must go to extreme values of these parameters, because in addition 
to needing to invert the natural hierarchy between yt and yb, one must overcome the big 
suppression for small fermion masses coming from the function T. 

Precisely the same argument applies to a broad class of Zee-type models, where the 
SM is supplemented by scalar multiplets whose weak isospin and hypercharge permit a 
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Yukawa coupling to the b quark and one of the other SM fermions. So long as the scalars 
do not mix and there are no new fermions to circulate in the loop, all such models have 
the same difficulty in explaining the Rb discrepancy. Below we will give some examples of 
models which, in contrast, are able to explain Ri,. 



5.3) Generalization to Nondiagonal Z Couplings 



We now turn to the more complicated case where mixing introduces off-diagonal cou- 
plings among the new particles. Because of mixing the couplings of the fermions to the Z 
will be matrices in the mass basis. Similar to eqs. (42) and (45) we write 



(85) 



where W^fn are the mixing matrices. An analogous expression gives the off-diagonal scalar- 
Z coupling in terms of the scalar mixing matrix, Us'^. Of course if all of the mixing 
particles share the same value for Is, then unitarity of the mixing matrices guarantees that 
the couplings retain this form in any basis. 

This modification of the neutral-current couplings has two important effects on the 
calculation of Sg'^. One is that the off-diagonal Z couplings introduce the additional graphs 
of the type shown in Fig. 5a,b, where the fermions or scalars on either side of the Z vertex 
have different masses. The other is that the mixing matrices spoil the relationship, eq. (81), 
whereby the term proportional to JF canceled in eq. (78). But this is only because of the 
mass-dependence of JF and A^. Therefore the cancellation still occurs if all of the particles 
that mix with each other are degenerate, as one would expect. Moreover the ultraviolet 
divergences still cancel since they are mass-independent. 

Evaluation of the graphs gives the following result at = M| = 0: 



S9l = 



327r2 



G 



diag 



(86) 



where Gdiag/327r^ represents the contribution involving only the diagonal Z couplings, and 
so is identical to the previously derived eq. (78). It is convenient to write it as 



G. 



diag 



nc \y 



f4> 



[2{g, - g^yfj^ir) + [-{g,y^ + gl + {gs)^^] (a^ - J^{r))]. (87) 



Here and in the following expressions we use the notation r = mj/M^ and r' = m^f,/M^. 
As before A<^ denotes the UV-divergent quantity A<^ = -|- 7 + ln(M|/47r//^) -|- ^. 
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The remaining terms in eq. (86) come from the new graphs of Fig. 5a,b, where the 
scalars or fermions on either side of the Z vertex have different masses, due to mixing: 



(88) 
(89) 



where J-L{r,r'), J-R{r,r') and Ts{x,x') are given by: 



^.(r-,r') 



1 



r ^ 1 / 
m r : m r 



r — 1 



r — 1 

= ix-l)\x'-l) + (x'-l)(x'-x) + ^V) + (x-l)(x-x'y 



Inr 



r' - 1 
r' - 1 



Inr' 



(90) 
(91) 



{x — l){x — x') 

(92) 



and x^ x' are the mass ratios x = M^/m'j and x' = M^,/m'j. These expressions have 
several sahent features which we now discuss. First, eqs. (87), (88) and (89) are obviously- 
much more complicated than eq. (82). In particular, it is no longer straightforward to 
simply read off the sign of the result. 

Second, the sum of the UV divergences in eqs. (87), (88) and (89), 



ff'4>4>' 



(93) 



is basis independent since a unitary transformation of the fields cancels between the Yukawa 
and neutral-current couplings. Thus it can be evaluated in the electroweak basis where the 
neutral-current couplings are diagonal and proportional to —g^ + 9l+ 9si which vanishes 
due to conservation of weak isospin and hypercharge at the scalar-fermion vertex. We are 
therefore free to choose the renormalization scale /i^ in ln(M|//x^) to take any convenient 
value. The M(^-dependence of A<^ makes G,^,^' look unsymmetric under the interchange of 
(j) and 0', but this is only an artifact of the way it is expressed. For example when there 
are only two scalars, G^,^/ is indeed symmetric under the interchange of their masses. 

Third, all the contributions except those of G<;6<^' ^^re suppressed by powers of ruf/M^ 
in the limit that the scalars are much heavier than the fermions. Thus to get a large 
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enough correction to g'^ requires that: (z) not all of the scalars be much heavier than the 
fermions which circulate in the loop, or (u) the scalars mix significantly and have the right 
charges so that G,^,^/ is nonnegligible and negative. We use option (u) in what follows to 
construct another mechanism for increasing Rb. 

Finally, even if the two fermions are degenerate, one does not generally recover the 
previous expression (78) that applied in the absence of mixing. This is because Dirac mass 
matrices are diagonalized by a similarity transformation, M — > uImUh, not a unitary 
transformation. The left- and right-handed mixing angles can differ even when the diag- 
onalized mass matrix is proportional to the identity. Thus, in contrast to eq. (93), the 
expression V f4,V*i i[{9l)^^ ~ (Qr)^^ ] is not invariant under transformations of the 

fields, because y^^ is rotated by Ur (recall that y^^ is the Yukawa coupling only for the 
RH /'s) whereas is rotated by Ul- 

We can get some insight into eqs. (88)-(92) by looking at special values of the pa- 
rameters. Let us assume there is a dominant Yukawa coupling y between the left-handed 
b quark and a single species of scalar and fermion, /i and 0i in the weak basis, 

>Cscalar = 2/01 /l7i-^ + h.C. (94) 

In the mass basis the couplings will therefore be 

y,, = yUlH^'Jr- (95) 

Now gauge invariance only relates the (1,1) elements of the neutral-current coupling ma- 
trices in the weak basis: 

{9s)''+9l-{9iy'^0. (96) 
There are three limiting cases in which the results become easier to interpret: 

• 1: If all the scalars arc degenerate with each other, and likewise for the fermions, then 
the nonmixing result of eq. (82) holds, except one must make the replacement 

gi-gi^ (Wh SmUl g, S^Ul - ^«)i\ (97) 

where Sm is the diagonal matrix of the signs of the fermion masses. 

• 2: If there are only two scalars and if they are much heavier than all of the fermions, 
only the term G<^<^/ is significant. Let 0i and 02 denote the weak-eigenstate scalars, and 
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and (j)' the mass eigenstates; then 



^9l-'^2^lt-lt)^'slFs{Ml/Ml)- (98) 
^s(r) = ^^^lnr-l, (99) 

where Cg and Sg are the cosine and sine of the scalar mixing angle. The function Fs{r) is 
positive except at r = 1 where it is zero, and so the sign of 5g\^ is completely controlled by 
the factor (/g^ — /g^). We see that to increase i?b it is necessary that Ig^ < /g^. 

• 3: When there are only two fermions, with weak eigenstates /i, /2 and mass eigenstates 
/, /' both much heavier than any of the scalars, then 



5g 



167r2 



11 2 , 22 2 11 



+ {9f - 9n)clslF^{m}/m},) - 2{gl^ - gl^) c.s.c^s^F.imj/m},)]; 



(100) 



where s^n and c^^ are the sine and cosine of the difference or sum of the LH and RH mixing 
angles, 9^ — Sm^R, depending on the relative sign Sm of the two fermion mass eigenvalues, 
and 

Fn{r) = Fs{r) = J. ^ Inr - 1; and F^r) = ^^Inr - 1. (101) 
l[r — 1) r — 1 

The function F^^ has some of the same properties as Fg = Fj^, including invariance under 
r — > 1/r, being positive semidefinite and vanishing at r = 1. Plots of these functions are 
shown in Figure 7. Note that the first line of eq. (100) is the same as (97). 

To get some idea of the error we have made by neglecting the mass of the Z boson one 
can compute the lowest order correction as in section 5.3. The answer is more complicated 
than for the case of diagonal Z couplings, except when the fermions are degenerate with 
each other and likewise for the bosons. In that case the answer is given again by eq. (84) 
except that gl {gl^^ and g[ {Un Sm ul g^ Sm Ul - gn^^, precisely as in eq. (97). 
Thus we would still expect it to be a small correction even when there is mixing of the 
particles in the loop. 

These simplifying assumptions can be used to gain a semi-analytic understanding 
of why certain regions of parameter space are favoured in complicated models, which is 
often missing in analyses that treat the results for the loop integrals as a black box. The 
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observations we make here may be useful when searching for modifications to a model 
that would help to explain i?^. The next two sections exemplify this by creating some 
new models that take advantage of our insights, and by elucidating previous findings in an 
already existing model, supersymmetry. 

5.4) Examples of Models That Work 

Besides ruling out certain classes of models, our general considerations also suggest 
what is required in order to explain R^. Obviously new fermions and scalars are required, 
whose Yukawa couplings allow them to circulate inside the loop. We give two examples, 
one with diagonal and one with nondiagonal couplings of the new particles to the Z boson. 

For our first example we introduce several exotic quarks F, P and N, and a new Higgs 
doublet 0, whose quantum numbers are listed in Table 5. The unorthodox electric charge 
assignments do not ensure cancellation of electroweak anomalies, but this can be fixed by 
adding additional fermions, like mirrors of those given, which do not contribute to Ri. 



Field 



Spin 



U^{1) 





Nr. 



1 

3 
3 

3 



2 
2 
1 

1 



q + l 

1 



Table 5 

Field Content and Charge Assignments: Electroweak quantum numbers for the new fields which are added 
to the SM to produce the observed value for Rf,. 

The hypercharges in Table 5 allow the following Yukawa interactions: 

jCy = y NuQ'A'e^j + gp PnFlH^e,, + (/n NuFlH'e,, + h.c, (102) 

where Cij is the 2x2 antisymmetric tensor, H is the usual SM Higgs doublet and = 
^ is the SM doublet of third generation LH quarks. When H gets its VEV, {H) = u, 
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we find two fermion mass eigenstates, p and n, whose masses are rup = QpV and = QnV 
and whose electric charges are Qp = q+1 and Qn = q- There are also two new scalar mass 
eigenstates, (p±, whose electric charges are Q+ = q + ^ and Q_ = q — ^. 

In the mass eigenstate basis, the Yukawa interactions with the new scalars are 

jCy = yniibL(p+ - ynntLip- + h.c, (103) 

from which we see that the n couples to the 6-quark as in eq. (76). 

The weak isospin assignments of the n are /^^ = ~ and I^^ = 0, so that g^ — g^ = —\- 
Therefore, from eq. (82), one obtains dg\^ < 0. The central value of Rb can be reproduced 
if 5g\ = —0.0067, which is easily obtained by taking y ~ 1 and r ^ 1, so that ^(r) ~ 1. 
The Yukawa coupling could be made smaller by putting the new scalars in a higher colour 
representation like the adjoint. 

We have not explored the detailed phenomenology of this model, but it is clearly not 
ruled out since wc arc free to make the new fermions and scalars as heavy as we wish. 
And since we can always take rUp = nin, there is no contribution to the oblique parameter 
T. The contribution to Rb does not vanish even as the masses become infinite, but this is 
consistent with decoupling in the same way as a heavy t quark, since the new fermions get 
their masses through electroweak symmetry breaking. The price we have to pay for such 
large masses is correspondingly large coupling constants. 

Next we build a model that uses our results for nondiagonal couplings to the Z. It is a 
simple modification of the SM that goes in the right direction for fixing the Rb discrepancy 
but not quite far enough in magnitude. Variations on the same theme can completely 
explain Rb at the cost of making the model somewhat more baroque. 

Our starting point is a two-Higgs doublet extension of the SM. We take the two Higgs 

fields, Hd = ( I and Hy_ = ( ) , to transform in the usual way under the SM gauge 

symmetry. It was explained earlier why this model does not by itself produce the desired 
effect, but eq. (98) suggests how to fix this problem by introducing a third scalar doublet, 

A = I ^_|_ 1 , which mixes with the other Higgs fields. The charge assignments of these 

fields, listed in Table 6, ensure that the two fields iJ+ and A+ can mix even though they 
have different eigenvalues for 73. 

In this model the new scalar field cannot have any Yukawa couplings to ordinary 
quarks since these are forbidden by hypercharge conservation. The only Yukawa couplings 
involving the LH 6-quark are those which also generate the mass of the t-quark: 



■"yuk — Ut 



V2t^^bH+ + h.c., (104) 
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Field 


Spin 


SUc{3) 






Hd 





1 


2 


1 

2 


Hu 





1 


2 


-4 


A 





1 


2 


+1 



Table 6 

Field Content and Charge Assignments: Electroweak quantum numbers for all of the scalars — including 
the SM Higgs doublet — of the three-doublet model. 

where yt = rrit/vu is the conventionally- normahzed Yukawa coupling. We imagine Vu to 
be of the same order as the single-Higgs SM value, and so we expect yt to be comparable 
to its SM size. 

The scalar potential for such a model very naturally incorporates — A"*" mixing. 
Gauge invariance permits quartic scalar interactions of the form X{HlA){Hl^Hd) + h.c, 
which generate the desired off-diagonal terms: X{A~^H^*v*Vd + A~^H^v^*) + h.c. 

Since the weak isospin assignments are I^^ = — | and I^" — -|-|, the colour factor 
is ric = 1 and the relevant Yukawa coupling is y = ytV^^ we see that eq. (98) predicts the 
following contribution due to singly-charged Higgs loops: 

Sgl = -:^2clslFs{r) , (105) 

with r being the ratio of the scalar mass eigenstates, r = M^/M^,. Taking optimistic 
values for the parameters'"^ {6s — f, 2c|s| — ^, yt — 1 and M^j^/Mfj,/ — 10), we find 
Sgb = -0.0043, which is two thirds of what is required: [dgD^^p = -0.0067 ± 0.0021. 

In addition to the contribution of the singly-charged scalar loops, one should consider 
those of the other nonstandard scalar fields we introduced. Since all of the scalars that 
mix have the same eigenvalue for /a, their contribution is given by eq. (82), which is small 
if the scalars are much heavier than the light fermions. Then only the t-quark contribution 
is important. In this limit there are appreciable contributions only from the three charged 

13 

Note that the charged-scalar mixing in this model is suppressed if one of the scalar masses gets very 
large compared to the weak scale. 
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scalar fields, one of which is eaten by the physical W boson and so is incorporated into 
the SM t-quark calculation, and the other two of which we have just computed. 

So, for an admittedly special region of parameter space, this simple model considerably 
ameliorates the Rb discrepancy, reducing it to a la effect. It is easy to adapt it so as to 
further increase 5g^ and also enlarge the allowed region of the model's parameter space. 
The simplest way is by increasing the size of the colour factor ric or the isospin difference 
Ig — /g . For instance the new scalar. A, could be put into a 4 of SUl{2) rather than a 
doublet, and be given weak hypercharge y = +|. Then the singly-charged state A+ has 
I^^ = — |, making /g — /g = —2, which is twice as big as for the doublet. More new 
scalars must be added to generate mixing amongst the singly-charged scalar states. 

A second variation would be let the two new Higgs doublets be colour octets since 
this gives more than a five-fold enhancement of 5g^ due to the colour factor ric = ^■ 
It is still possible to write down quartic scalar interactions which generate the desired 
scalar mixings. Either of these models has much more room to relax the previously tight 
requirements for optimal scalar masses and mixings. 

5.5) The Supersymmetric Case 

Let us now apply the above results to gain some insight into what would be necessary 
to explain Ri, in supersymmetric extensions of the standard model. There are two kinds 
of contributions involving the top-quark Yukawa coupling, which one expects to give the 
dominant effect. These are the couplings of the left-handed b quark to the second Higgs 
doublet and the top quark, or to the corresponding Higgsinos and top squarks. 



Of these, the second one gives a loop contribution like that of the two-Higgs doublet 
models discussed above: it has the wrong sign for explaining R^. Since the mass of the 
charged Higgs is a free parameter in supersymmetric models, we can imagine making it 
large enough compared to nit so that, according to eq. (82), it has only a small effect on 
Rt,. We therefore concentrate on the Higgsino-squark part. The charged Higgsino mixes 
with the Wino, and the right-handed top squark mixes with its chiral counterpart, so in 
the notation of (94), we have /i = /12 7 /2 = W, (f)i — and 02 — i^^- The corresponding 
charge matrices for the couplings to the are 



(106) 




(107) 
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Because there are two possible colour combinations for the internal lines of the loops 
diagram, the colour factor in eqs. (87)- (89) is ric = 2. 

Before exploring the full expression for Sg^ we can discover what parameter ranges 
are the most promising by looking at the limiting cases described by eqs. (97)-(100). 
The most important lessons from these approximations follow from the charge matrices 
(107). We do not want the squarks to be much heavier than the charginos because then 
eq. (98) would apply and give the wrong sign for the correction due to the sign of the 
isospin difference between the squarks. The other two cases, where the squarks are not 
much heavier than the charginos, manifest a strong suppression of the result unless the 
chargino mixing angles are such that sin(^j^ — SmOa) is large, where Sm is the sign of the 
determinant of the chargino mass matrix. If on the other hand sin(^^ — Sm^ji) = 0, there is 
exact cancellation between g^^ and g^ in these equations because of the fact that g^^ = g^, 
for the charginos. In summary, our analytic formulas indicate that the favoured regions of 
parameter space for increasing Ri, are where 

tan^j^ tan^j^ = —Sm = — sign(m/ m//), (108) 

and at least one of the squarks is not much heavier than the charginos. 

In supersymmetric models the Yukawa coupling that controls the largest contribution 
to is that of the top quark, and it depends on the ratio of the two Higgs VEV's, 
tan j3 = V2/V1, by 

where v = {vj + ^2)"^^^ = 174 GeV. Therefore it is important to find tan/? in terms of the 
chargino masses and mixing angles. The chargino mass matrix is given by 

A* 9V2 \ f rrif ^ \i( = ( ^LCRmf + SLSnrufr SnCLUif - 0^8^171 f 

gVi M2 ) \ rUf ) \CnSLmf - SnCLrUf SLSnTUf + CLCnrUf 

(110) 

where ^ is the coefficient of H1H2 in the superpotential and M2 is the soft-SUSY-breaking 
mass term for the Wino. It follows that 

.^^^^ruftane -mf.tane (111) 
m/ tan^i — uif tandn 

The above considerations allow us to understand why values of tan /3 near unity are 
necessary for a supersymmetric solution to the Rb problem. From eq. (Ill) and the 
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maximization condition (108) we see that tan/3 is restricted to lie between |m//m//| and 
\mfi /mf\. Eq. (108) together with (110) also implies 

\ = \/2Mw sin (3; cl\ \ = V2M^cosP, (112) 

This means that average value of the two chargino masses can be no greater than M^y, so 
that the ratio |m//m/'| cannot differ much from unity unless one of the charginos is much 
lighter than the W boson. Using the LEP 1.5 limit of 65 GeV for the lightest chargino [26] 
this would then require that tan/? < 1.5. 




ure 8: The dependence of (7^ on the various supersymmetric parameters. Since (7^ depends only on mass 
ratios in our approximation, the units of mass are arbitrary, with the masses of all the charginos and 
squarks which are not being varied set to unity. 

In the case that none of our simplifying limits apply, we have searched the parameter 
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space of the three independent ratios between the two scalar masses and the two fermion 
masses, and the three mixing angles 9^^ 6^, 9s to find which regions are favourable for 
increasing i?^. Figures 8a-d show the shift in as a function of pairs of these parameters, 
using the Yukawa coupling (109) corresponding to a top quark mass of 174 GeV and 
the theoretical preference for tan/? > 1 (we implement the latter by setting = for 
parameters that would give tan/3 < 1). As shown in Table 1, one needs Sg'l = —0.0067 
in order to explain the observed value of The values of the masses are taken to be 
Mff) = Mff)' — ruf mf = 1 (in arbitrary units), except for those that are explicitly 
varied in each figure. In Fig. 8a we look at the situation in which tan^^ = tan^^. = 1, in 
contradiction to condition (108), and vary the scalar mixing angle and the mass of mostly- 
scalar in the limit of zero squark mixing. The sign of (7^ has the wrong value, as predicted 
by eq. (98). Fig. 8b shows the same situation except that now tan^j, = — tan^^ = 1, in 
accordance with eq. (108). Then the sign of g'^ is negative, as desired, and has the right size 
for substantial ranges of 9s and M^. In Fig. 8c we keep all the masses nearly degenerate 
and set = to show the dependence on tan^^ and tan^n- It is easy to see that g'l 
has the correct sign and largest magnitude (which is also almost as large as needed) when 
condition (108) is satisfied. Finally in Fig. 8d we show the dependence on the masses of 
the mostly- Wino fermion and on 9^ when = and tan^^ = —1, showing again the 
preference for mixing angles obeying (108), as well as some enhancement when there is a 
hierarchy between the two chargino masses. 



One might therefore get the impression that it is easy to explain R(, using super- 
symmetric contributions to the Zbb vertex. The problem is that to get a large enough 
contribution one is driven to a rather special region of parameter space, which comes 
close to satisfying condition (108). As mentioned above, the consequent condition (112) 
prevents one from making the chargino masses arbitrarily heavy. This, coupled with the 
suppression in i?& when the squarks are heavier than the charginos, means that all the 
relevant supersymmetric particles must be relatively light, except the charged Higgs which 
has to be heavy to suppress the wrong-sign contribution from H'^-t loops. Thus in the 
example of Fig. 8c, the preferred values of = 1, s^, = ±1, = c^, = imply that 
nif = V sin P and ruf = vcos/3, while = M2 = 0, which are precisely the circumstances 
of the supersymmetric models considered in Refs. [27] and [28] . Fig. 8d, on the other hand, 
has its maximum value of Rb at 0^ = 811 = 0^ = —s^ = 1, implying tan/3 = 1 and thus 
from (112) that |m//| -|- |m/| = 2My^^. Because the lightest chargino mass is constrained 
by experimental lower limits, there is little parameter space for getting a large hierarchy 
between the two chargino masses, as one would want in the present example in order to 
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get the full shift of —0.0067 in g^. Our analysis allows one to pinpoint just where the 
favorable regions are for solving the Ri, problem. 

We thus see that it is possible to understand many of the conclusions in the literature 
[27]-[31] on supersymmetry and Rb using some rather simple analytic formulas. These 
include the preference for small values of tan /? as well as light higgsinos and squarks. 



6. Future Tests 



If we exclude the possibility that the experimental value of Rb is simply a 3.7cr sta- 
tistical fluctuation, we can expect that, once the LEP collaborations have completed their 
analyses of all the data collected during the five years of running at the Z pole, the ^R^ 
crisis' will become an even more serious problem for the standard model. (Of course, it is 
wise to keep in mind that there may be a simple explanation, namely that some systematic 
uncertainties in the analysis of the experimental data are still not well understood or have 
been underestimated.) In sections 3-5 we have discussed a variety of models of new physics 
which could account for the experimental measurement of Rb- The next obvious step is to 
consider which other measurements may be used to reveal the presence of this new physics. 

The most direct method of finding the new physics is clearly the discovery of new 
particles with the correct couplings to the Z and the b quark. However, failing that, there 
are some indirect tests. For example, many of the new-physics mechanisms which have 
been analysed in this paper will affect the rate for some rare B decays in a predictable way. 
The rates for the rare decays B Xs£~^i~ and B XgUiy are essentially controlled by 
the Zbs effective vertex F'^^ , since additional contributions (such as box diagrams and Z-7 
interference) are largely subleading.^^ In the SM, in the approximation made throughout 
this paper of neglecting the 6-quark mass and momentum, a simple relation holds between 
the dominant rrit vertex effects in Rb and in the effective Zbs vertex F'^ : 



6. |^^^|2 



An additional constraint is that the Hghtest Higgs boson mass m^o vanishes at tree level when 

tan/3=l, and a very large splitting between the top squark masses is needed for the one-loop corrections 

to m^o to be large enough. This is why ref. [29] finds less than the desired shift in Rb in the minimal 

supersymmetric standard model. We thank J. Lopez for clarifying this point. 
15 

Due to the absence of Z-7 interference and of large renormalization-group-induced QCD corrections, 
the process B^X^ vv represents theoretically the cleanest proof of the effective Zhs vertex [32] . 



53 



where ST^^'^'^ is defined as in eq. (26) with the SM form factor as given in eqs. (27) and (38). 
The meaning of eq. (113) is that, within the SM, the Zbs effective vertex measurable in Z- 
mediated B decays represents a direct measurement of the mt-dependent vertex corrections 
contributing to i?5, modufo a ratio of the relevant CKM matrix elements. In particular, 
both corrections vanish in the rrit limit. The question is now: how is this relation 
affected by the new physics invoked in Sees. 3-5 to explain Rb7 

Consider first the tree-level b-b' mixing effects analysed in Sec. 3. It is straightforward 
to relate the corrections of the LH and RH Zbb couplings to new tree-level mixing-induced 
FCNC couplings gl\. In this case eq. (5) reads 

= Y.^C).,nK^nK%M'^' . (114) 
aw 

Hence g^j^j^ involve the same gauge couplings and mixing matrices that determine the 
deviation from the SM of the flavour-diagonal b couplings. 

It is also true that, for many models of new physics, the loop corrections to the Zbb 
vertex would change the effective Zbs vertex in much the same way, therefore inducing 
computable modifications to the SM electroweak penguin diagrams. In these models, for 
each loop diagram involving the new states /, /' and their coupling to the 6-quark gffh-, 
there will be a similar diagram contributing to F^^ that can be obtained by the simple 
replacement gffh 9ff's- For example, the general analysis of t-quark mixing effects 
presented in Sec. 4 can be straightforwardly applied to Z-mediated B decays. Deviations 
from the SM predictions for the B Xs£+£~ and B ^ X sVi' decay rates can be easily 
evaluated by means of a few simple replacements like |Vt6p — > V^^Vts and |Vt'6p — > V^^Vfa 
in all our equations. To a large extent, this is also true for SUSY models. Indeed, the 
analysis of the SUSY contributions to the Zbs form factor [34] can teach much about SUSY 
effects in i?^. And once a particular region of parameter space suitable to explain the i?& 
problem is chosen, a definite numerical prediction for the B Xs£~^£~ and B Xsvy 
decay rates can be made. 

This brief discussion shows that, for a large class of new-physics models, the new 
contributions to i?6 and to the effective V'^^ vertex are computable in terms of the same 
set of new-physics parameters. Therefore, for all these models, the assumption that some 
new physics is responsible for the deviations of Rh from the SM prediction will imply a 
quantitative prediction of the corresponding deviations for Z-mediated B decays. 

For example, the particular case of mixing of the top-quark with a new isosinglet T', and the corre- 
sponding effects induced on the Zhs vertex, was studied in Ref. [33] through an analysis very similar to 
that of Sec. 4. 
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However, this statement cannot be applied to all new-physics possibilities. For ex- 
ample, if a new Z' boson is responsible for the measured value of then no signal can 
be expected in B decays, since in this case the new physics respects the GIM mechanism. 
This would also be true if mfo-dependent effects are responsible for the observed deviations 
in i?5 as could happen, for example, in the very large tan /? region of multi-Higgs-doublet 
or SUSY models. More generally, the loop contributions of the new states /, /' can be 
different, since gffs is not necessarily related to g^fh-, and in particular, whenever the 
new physics involved in couples principally to the third generation, it is quite possible 
that no sizeable effect will show up in B decays. Still, the study of S ^ Xsl'^l~ and 
B ^ XgVv could help to distinguish between models that do or do not significantly affect 
these decays. 

Unfortunately, at present only upper limits have been set on the branching ratios for 
B — > Xsl'^l~ [35]-[37] and B — > XgVv [32]. Since these limits are a few times larger than 
the SM predicitons, they cannot help to pin down the correct solution to the Kh problem. 
However, future measurements of these rare decays at B factories could well confirm that 
new physics is affecting the rate of 6-quark production in Z decays, as well as give some 
hints as to its identity. If no significant deviations from the SM expectations are detected, 
this would also help to restrict the remaining possibilities. 

7. Conclusions 

Until recently, the SM has enjoyed enormous success in explaining all electroweak 
phenomena. However, a number of chinks have started to appear in its armour. There are 
currently several disagreements between theory and experiment at the 2cr level or greater. 
They are: Rh = Ffe/Fhad (3.7cr), = Fc/Fhad (2.5cr), the inconsistency between as 
measured at LEP with that determined at SLC (2.4a), and ^°b(t) (2. Oct). Taken together, 
the data now exclude the SM at the 98.8% confidence level. 

Of the above discrepancies, it is essentially only i?6 which causes problems. If Rh by 
itself is assumed to be accounted for by new physics, then the fit to the data despite the 
other discrepancies is reasonable (x^^^/d.o.f. = 15.5/11) - the other measurements could 
thus be regarded simply as statistical fluctuations. 

In this paper we have performed a systematic survey of new-physics models in order 
to determine which features give corrections to Rh of the right sign and magnitude. The 
models considered can be separated into two broad classes: those in which new Zhh cou- 
plings appear at tree level, by Z or 6-quark mixing with new particles, and those which give 
loop corrections to the Zhh vertex. The latter type includes t-quark mixing and models 
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with new scalar s and fermions. We did not consider technicolour models or new gauge 
bosons appearing in loops since these cases are much more model-dependent. 

The new physics can modify either the left-handed or right-handed Zbb couplings, 
or g^. To increase Rb to its experimental value, must be negative and have a 
magnitude typical of a loop correction with large Yukawa couplings. Thus 5g^ could either 
be a small tree-level effect, or a large one-loop effect. On the other hand, the SM value of 
g^ is opposite in sign to its LH counterpart and is about five times smaller. Therefore one 
would need a large tree-level modification to g^ to explain for i?^. 

Here are our results: 

• (1) Tree-level Effects: It is straightforward to explain if the Z or b mix with new 
particles. With Z-Z' mixing there are constraints from neutral-current measurements, but 
these do not exclude all models. Using b-b' mixing is easier since the experimental value 
of Rh can be accommodated by bj^-b'j^ or b^-b'^^ mixing. If the mixing is in the LH b sector, 
then solutions are possible so long as < —1/2. An additional possibility with /g^ > 
and very large LH mixing, though perhaps unappealing, is still viable. For RH b mixing, 
if /g^ > then small mixing is permitted, while if I^^ < 0, large mixing is necessary. 
Interestingly, the required large 6-mixing angles are still not ruled out phenomenologically. 
A number of papers in the literature have appealed to b-b' mixing to explain Ri,. Our 
"master formula" (8) and Table 2 include all of these models, as well as many others. 

• (2) Loops: t-t' Mixing: In the presence of t-t' mixing, the SM radiative correction can 
be reduced, depending on the weak isospin quantum numbers of the t' as well as on the 
LH and RH mixing angles. However, we found that it is not possible to completely explain 
Rb via this method. The best we can do is to decrease the discrepancy between theory 
and experiment to about 2a. Such a scenario predicts the existence of a light (~ 100 GeV) 
charge 2/3 quark, decaying primarily to Wb. 

• (3) Loops: Diagonal Couplings to the Z: We considered models with exotic fermions and 
scalars coupling to both the Z and 6-quark. We assumed that the couplings to the Z are 
diagonal, i.e. there are no flavour-changing neutral currents (FCNC's). The correction 5g^^ 
can then be written in a simple form, eq. (82). The key point is that 5g\^ is proportional to 
ll^ — iIj^, where ^ ^ is the third component of weak isospin of the fermion field fi^^j^ in 
the loop. This explains at a glance why many models, such as multi-Higgs-doublet models 
and Zee- type models, have difficulty explaining Rf,. Since the dominant contributions in 
these models typically have top-type quarks (/g^ = |, Ig^ = 0) circulating in the loop, 
they give corrections of the wrong sign to Ri,. However, these considerations did permit 
us to construct viable models of this type which do explain Rf,. Two such examples are 
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given Sec. 5.4, and many others can be invented. 

• (4) Loops: Nondiagonal Couplings to the Z : We also examined models with exotic 
fermions and scalars which were allowed to have nondiagonal couplings to the Z. Such 
FCNC's can occur when particles of different weak isospin mix. The correction 6g\^ is 
much more complicated (eq. (86)) than in the previous case; even its sign is not obvious. 
However there are several interesting limiting cases where it again becomes transparent. 
The contributions to Rh of supersymmetry fall into this category, which we discussed in 
some detail. 
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Note Added: After completing this work we became aware of ref. [38], which dis- 
cusses a different region of parameter space in SUSY models than the one we focused on. 
Because of our criterion of explaining the entire Rb discrepancy rather than only reducing 
its statistical significance, we exclude the region in question. 
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